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Abstract 

We propose a cosmological scenario in which the hot big bang universe is 
produced by the collision of a brane in the bulk space with a bounding orb- 
ifold plane, beginning from an otherwise cold, vacuous, static universe. The 
model addresses the cosmological horizon, flatness and monopole problems 
and generates a nearly scale-invariant spectrum of density perturbations with- 
out invoking superluminal expansion (inflation). The scenario relies, instead, 
on physical phenomena that arise naturally in theories based on extra di- 
mensions and branes. As an example, we present our scenario predominantly 
within the context of heterotic M-theory. A prediction that distinguishes this 
scenario from standard inflationary cosmology is a strongly blue gravitational 
wave spectrum, which has consequences for microwave background polariza- 
tion experiments and gravitational wave detectors. 
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I. INTRODUCTION 



The Big Bang model provides an accurate account of the evolution of our universe from 
the time of nucleosynthesis until the present, but does not address the key theoretical puz- 
zles regarding the structure and make-up of the Universe, including: the flatness puzzle 
(why is the observable universe so close to being spatially flat?); the homogeneity puzzle 
(why are causally disconnected regions of the universe so similar?); the inhomogeneity puz- 
zle (what is the origin of the density perturbations responsible for the cosmic microwave 
background anisotropy and large-scale structure formation? And why is their spectrum 
nearly scale-invariant?); and the monopole problem (why are topological defects from early 
phase transitions not observed?). Until now, the leading theory for resolving these puzzles 
has been the inflationary model of the universe.BH The central assumption of any infla- 
tionary model is that the universe underwent a period of superluminal expansion early in 
its history before settling into a radiation-dominated evolution. Inflation is a remarkably 
successful theory. But in spite of twenty years of endeavor there is no convincing link with 
theories of quantum gravity such as M-theory. 

In this paper, we present a cosmological scenario which addresses the above puzzles but 
which does not involve inflation. Instead, we invoke new physical phenomena that arise 
naturally in theories based on extra dimensions and branes. Known as "brane universe" 
scenarios, these ideas first appeared in Refs. || and 0. However, only recently were they given 
compelling motivation in the work of Hof ava and Witteni and in the subsequent construction 
of heterotic M-theory by Lukas, Ovrut and Waldram.i Complementary motivation was 
provided both in superstring theory&H and in non-string contexts.0'0 Many of the ideas 
discussed here are applicable, in principle, to any brane universe theory. For example, in 
discussing the features of our model, we draw examples both from the Randall- Sundrum 
mo del0 and from heterotic M-theory. However, in this paper, we emphasize heterotic M- 
theory. This is done for specificity and because, by doing so, we know that we are working 
in a theory that contains all the particles and interactions of the standard model of particle 
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physics. Hence, we are proposing a potentially realistic theory of cosmology. 

Specifically, our scenario assumes a universe consisting of a five- dimensional space-time 
with two bounding (3+l)-dimensional surfaces (3-branes) separated by a finite gap spanning 
an intervening bulk volume. One of the boundary 3-branes (the "visible brane") corresponds 
to the observed four-dimensional universe in which ordinary particles and radiation propa- 
gate, and the other is a "hidden brane." The universe begins as a cold, empty, nearly BPS 
(Bogolmon'yi-Prasad-Sommerfieldlii) ground state of heterotic M-theory, as described by 
Lukas, Ovrut and Waldram.i The BPS property is required in order to have a low-energy 
four- dimensional effective action with M = 1 supersymmetry. The visible and hidden branes 
are flat (Minkowskian) but the bulk is warped along the fifth dimension. 

In addition to the visible and hidden branes, the bulk volume contains an additional 
3-brane which is free to move across the bulk. The bulk brane may exist initially as a 
BPS state, or it may spontaneously appear in the vicinity of the hidden brane through a 
process akin to bubble nucleation. The BPS condition in the first case or the minimization 
of the action in the second case require that the bulk brane be flat, oriented parallel to the 
boundary branes, and initially at rest. Non-perturbative effects result in a potential which 
attracts the bulk brane towards the visible brane. We shall assume that the bulk brane is 
much lighter than the bounding branes, so that its backreaction is a small correction to the 
geometry. See Figure 1. 

The defining moment is the creation of the hot big bang universe by the collision of the 
slowly moving bulk brane with our visible brane. Although the universe may exist for an 
indefinite period prior to the collision, cosmic time as normally defined begins at impact. The 
bulk and visible branes fuse through a "small instanton" transition, during which a fraction 
of the kinetic energy of the bulk brane is converted into a hot, thermal bath of radiation and 
matter on the visible brane. The universe enters the hot big bang or Friedmann-Robertson- 
Walker (FRW) phase. Notably, instead of starting from a cosmic singularity with infinite 
temperature, as in conventional big bang cosmology, the hot, expanding universe in our 
scenario starts its cosmic evolution at a finite temperature. We refer to our proposal as 
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the "ekpyrotic universe", a term drawn from the Stoic model of cosmic evolution in which 
the universe is consumed by fire at regular intervals and reconstituted out of this fire, a 
conflagration called ekpyrosis.0 Here, the universe as we know it is made (and, perhaps, 
has been remade) through a conflagration ignited by collisions between branes along a hidden 
fifth dimension. 

Forming the hot big bang universe from colliding branes affects each of the cosmological 
problems of the standard big bang model. First, the causal structure of space-time in the 
scenario differs from the conventional big bang picture. In standard big bang cosmology, 
two events separated by more than a few Hubble radii are causally disconnected. This 
relationship between the Hubble radius and the causal horizon applies to FRW cosmologies 
that are expanding subluminally, as in the standard hot big bang model, but does not apply 
to more general cosmologies, such as de Sitter space-time or the scenario we will describe. 
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FIG. 1. One possible set of initial conditions of the ekpyrotic scenario has the universe beginning 
in a cold, vacuous, nearly BPS state consisting of two static massive orbifold planes and a warped 
geometry in the intervening bulk (a) in which the curvature is low near the rightmost orbifold plane 
(the hidden brane) and high near the leftmost orbifold plane (the visible brane). Spontaneously, a 
bulk brane peels away from the hidden brane over some region of space (b) , forming a terrace. The 
edges of the terrace expand outwards at light speed, while the interior moves very slowly towards 
the opposing visible brane. Although the bulk brane is flat on average, quantum fluctuations 
produce ripples over a wide range of length scales as the brane traverses the bulk (c). When 
the bulk brane collides with the visible brane, the ripples result in different regions colliding and 
reheating at slightly different times (d), thereby impressing a spectrum of density fluctuations on 
the visible universe. The energy from the collision is translated into matter and radiation, heating 
the universe to a temperature a few orders of magnitude smaller than the unification scale. 

In our scenario, the collision sets the initial temperature and, consequently, the Hubble 
radius at the beginning of the FRW phase. The Hubble radius is generally infinitesimal 
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compared to the collision region. Two events outside the Hubble radius are correlated since 
local conditions have a common causal link, namely, the collision with the bulk brane. Here 
we take advantage of the fact that the brane is a macroscopic, non-local object and exists 
for an indefinitely long period prior to the collision. That is, there is no direct connection 
between the time transpired preceding the collision or causality and the Hubble time. This 
feature provides a natural means for resolving the horizon problem. 

As we have noted above, the boundary branes and the bulk brane are initially flat 
and parallel, as demanded by the BPS condition. Furthermore, the motion of the bulk 
brane along the fifth dimension maintains flatness (modulo small fluctuations around the 
flat background). Hence, the hot big bang universe resulting from the collision of a flat 
bulk brane with a flat visible brane is spatially flat. In other words, we address the flatness 
problem by beginning near a BPS ground state. 

We do not require the initial state to be precisely BPS to resolve the horizon and flatness 
problems. It suffices if the universe is flat and homogeneous on scales ranging up to the 
(causal) particle horizon, as should occur naturally beginning from more general initial 
conditions. In the ekpyrotic scenario, the distance that particles can travel before collision is 
exponentially long because the bulk brane motion is extremely slow. As a result, the particle 
horizon at collision can be many more than 60 e-folds larger than the Hubble radius at 
collision (where the latter is determined by the radiation temperature at collision). That is, 
rather than introducing superluminal expansion to resolve the horizon and flatness problems, 
the ekpyrotic model relies on the assumption that the universe began in an empty, quasi- 
static BPS state which lasted an exponentially long time prior to the beginning of the hot 
big bang phase. 

Quantum fluctuations introduce ripples in the bulk brane as it moves across the fifth 
dimension. During this motion there is a scale above which modes are frozen in, and below 
which they oscillate. This scale decreases with time in a manner akin to the Hubble radius 
of a collapsing Universe. The fluctuations span all scales up to this freeze-out scale, and we 
assume they begin in their quantum mechanical ground state. As the brane moves across the 
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bulk, the effective Hubble radius shrinks by an exponential factor while the wavelengths of 
the modes decrease only logarithmically in time. Consequently, modes that begin inside the 
initial freeze-out scale end up exponentially far outside it, and exponentially far outside the 
final Hubble radius at the time of collision. The ripples in the bulk brane cause the collision 
between the bulk brane and the visible brane to occur at slightly different times in different 
regions of space. The time differences mean regions heat up and begin to cool at different 
times, resulting in adiabatic temperature and density fluctuations. Hence, as interpreted by 
an observer in the hot big bang FRW phase, the universe begins with a spectrum of density 
fluctuations that extend to exponentially large, super-horizon scales. 

Of course, the spectrum of energy density perturbations must be nearly scale-invariant 
(Harrison-Zel'dovichB) to match observations of large-scale structure formation and tem- 
perature anisotropics in the cosmic microwave background. We find that this condition is 
satisfied if the potential attracting the bulk brane towards the visible brane is ultra weak at 
large separations. This is consistent with potentials generated by non-perturbative effects 
such as the exchange of virtual M2-branes (wrapped on holomorphic curves) between the 
bulk brane and either of the boundary branes. 

While the spectrum of perturbations is approximately scale invariant, as with inflation 
there are small deviations from scale invariance. In the examples considered here, the 
spectrum is blue (the amplitude increases as the wavelength decreases), in contrast to typical 
inflationary models. With exponentially flat potentials, the spectrum is only marginally 
blue, consistent with current observations. On the other hand, the potential has no effect 
on the tensor (gravitational wave) perturbations, so the tensor spectrum is strongly blue 
(spectral index tit ~ 2), in contrast to the slightly red (jit < 0) spectrum predicted in most 
inflationary models. This prediction may be tested in near-future microwave background 
anisotropy and gravitational wave detector experiments. 

For some aspects of the ekpyrotic scenario, such as the generation of quantum fluctua- 
tions, the description from the point-of-view of an observer on the bulk brane is the most 
intuitive. For that observer, the scale factor and Hubble radius appear to be shrinking 
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because the warp factor decreases as the bulk brane moves across the fifth dimension. How- 
ever, as observed from the near-stationary boundary orbifold planes, the universe is slowly 
expanding due to the gravitational backreaction caused by the bulk brane motion. Indeed, 
a feature of the scenario is that the bulk brane is responsible for initiating the expansion 
of the boundary branes. Furthermore, as we shall show, the brane gains kinetic energy due 
to its coupling to moduli fields. Upon impact, the bulk brane is absorbed by the visible 
brane in a so-called small-instanton phase transition (see Sec. II for details). This transi- 
tion can change the gauge group on the visible brane. (For example, before collision the 
gauge group on the visible brane might be one of high symmetry, such as E 6 , whereas after 
collision it becomes the standard model gauge group.) Furthermore, the number of light 
families of quarks and leptons on the visible brane may change during the transition. (For 
example, the visible brane might make a transition from having no light families of quarks 
and leptons to having three.) Upon collision, the kinetic energy gained by the bulk brane 
is converted to thermal excitations of the light degrees of freedom, and the hot big bang 
phase begins. Hence, the brane collision is not only responsible for initiating the expansion 
of the universe, but also for spontaneously breaking symmetries and for producing all of the 
quarks and leptons. 

If the maximal temperature lies well below the mass scale of magnetic monopoles (and 
any other cosmologically dangerous massive, stable particles or defects), none will be gen- 
erated during the collision and the monopole problem is avoided. 

Although the ideas presented in this work may be applicable to more general brane- 
world scenarios, such as Randall-Sundrum, in developing our scenario we have felt that it 
is important to take as a guiding principle that any concepts introduced in this scenario 
be consistent with string theory and M-theory. By founding the model on concepts from 
heterotic M-theory, one knows from the outset that the theory is rich enough to contain 
the particles and symmetries necessary to explain the real universe and that nothing we 
introduce interferes with a fundamental theory of quantum gravity. We emphasize that our 
scenario does not rely on exponential warp factors (which are inconsistent with BPS ground 
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states in heterotic M-theory) nor does it require large (millimeter-size) extra dimensions. For 
example, we consider here a bulk space whose size is only four or five orders of magnitude 
larger than the Planck length, consistent with Hof ava-Witten phenomenology .EBB All brane 
universe theories, including heterotic M-theory, suffer from some poorly understood aspects. 
For example, we have nothing to add here about the stability of the final, late time vacuum 
brane configuration. We will simply assume that branes in the early universe move under 
their respective forces until, after the big bang, some yet unknown physics stabilizes the 
vacuum. 

Figure 1 summarizes the conceptual picture for one possible set of initial conditions. The 
remainder of this paper discusses our attempt to transform the conceptual framework into 
a concrete model. For this purpose, a number of technical advances have been required: 

• an understanding of the perturbative BPS ground state (Section II) and how it can 
lead to the initial conditions desired for our scenario (Section III); 

• a moduli space formulation of brane cosmology (Section IV A); 

• a derivation of the equations of motion describing the propagation of bulk branes in a 
warped background in heterotic M-theory, including non-perturbative effects (Section 
IV B-D); 

• a computation of the bulk brane-visible brane collision energy, which sets the initial 
temperature and expansion rate of the FRW phase (Section IV D); 

• an analysis of how the gravitational backreaction due to the motion of the bulk brane 
induces the initial expansion of the universe (Section IV E); 

• a computation of how ripples in the bulk brane translate into density perturbations 
after collision with the visible brane (Section V A); 

• a theory of how quantum fluctuations produce a spectrum of ripples on the bulk brane 
as it propagates through a warped background (Section V B); 
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• a determination of the generic conditions for obtaining a nearly scale invariant spec- 
trum and application of general principles to designing specific models (Section V 
B); 

• a calculation of the tensor (gravitational wave) perturbation spectrum (Section V C); 

• a recapitulation of the full scenario explaining how the different components rely on 
properties of moduli in 5d (Section VI A); 

• a fully worked example which satisfies all cosmological constraints (Section VI A); 

• and, a comparison of the ekpyrotic scenario with inflationary cosmology, especially 
differences in their predictions for the fluctuation spectrum (Section VI B). 

In a subsequent paper, we shall elaborate on the moduli space formulation of brane cos- 
mology and show how it leads to a novel resolution of the singularity problem of big bang 
cosmology.0 

The ekpyrotic proposal bears some relation to the pre-big bang scenario of Veneziano et 
a/.&i^l which begins with an almost empty but unstable vacuum state of string theory but 
which, then, undergoes superluminal deflation. Several important conceptual differences are 
discussed in Section VIB. Models with brane interactions that drive inflation followed by 
brane collision have also been considered.0~il Applications of the moduli space of M-theory 
and Hofava-Witten theory to cosmology have been explored previously in the context of 
inflation&il The distinguishing feature of the ekpyrotic model is that it avoids inflation 
or deflation altogether. A non-inflationary solution to the horizon problem was suggested 
in Ref. [29], but it is not clear how to generate a nearly scale-invariant spectrum of density 
fluctuations without invoking inflation. 

We hope that our technical advances may be useful in exploring other variants of this 
scenario. Some aspects remain more speculative, especially the theory of initial conditions 
(as one might expect) and non-perturbative contributions. A detailed understanding of 
these latter aspects awaits progress in heterotic M-theory. 
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II. TERMINOLOGY AND MOTIVATION FROM HETEROTIC M-THEORY 

In this section, we briefly recount key features of heterotic M-theory that underlie and 
motivate the example of our ekpyrotic scenario given in this paper. Those who wish to 
understand the basic cosmological scenario without regard to the heterotic M-theory context 
can proceed directly to the next section. 

Heterotic M-theory has its roots in the work of Hof ava and Wittenl who showed that com- 
pactifying (eleven-dimensional) M-theory on an S1/Z2 orbifold corresponds to the strong- 
coupling limit of heterotic E s x E s (ten-dimensional) string theory. Compactifying an ad- 
ditional six dimensions on a Calabi-Yau three- fold leads, in the low energy limit, to a four- 
dimensional Af = 1 supersymmetric theory ,0 the effective field theory that underlies many 
super symmetric theories of particle phenomenology. By equating the effective gravitational 
and grand-unified coupling constants to their physical values, it was realized that the orbifold 
dimension is a few orders of magnitude larger than the characteristic size of the Calabi-Yau 
space. EE There is, therefore, a substantial energy range over which the universe is effec- 
tively five-dimensional, being bounded in the fifth-dimension by two (3 + l)-dimensional 
"end-of-the- world" S /Z2 orbifold fixed planes. By compactifying Hofava-Witten theory 
on a Calabi-Yau three-fold, in the presence of a non-vanishing G-flux background (that 
is, a non-zero field strength of the three-form of eleven dimensional M-theory) required by 
anomaly cancellation, the authors of Refs. |6] and ^ were able to derive the explicit effective 
action describing this five-dimensional regime. This action is a specific gauged version of 
M = 1 supergravity in five dimensions and includes "cosmological" potential terms that al- 
ways arise in the gauged context. It was shown in Ref. |6] that these potentials support BPS 
3-brane solutions of the equations of motion, the minimal vacuum consisting of two 3-branes, 
each coinciding with one of the S1/Z2 orbifold fixed planes. These boundary 3-branes (the 
visible brane and the hidden brane) each inherit a (spontaneously broken) Af = 1 Eg su- 
per gauge multiplet from Hofava-Witten theory. This five-dimensional effective theory with 
BPS 3-brane vacua is called heterotic M-theory. It is a fundamental paradigm for "brane 
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universe" scenarios of particle physics. 

In order to support a realistic theory, heterotic M-theory must include sufficient gauge 
symmetry and particle content on the 3-brane boundaries. In the compactification discussed 
above, the authors of Refs. |6| and |3(] initially made use of the standard embedding of the spin 
connection into the gauge connection, leading to an E 6 gauge theory on the visible 3-brane 
while the gauge theory on the hidden 3-brane is Eg. Subsequently, more general embeddings 
than the standard one were considered.0 Generically, such non-standard embeddings (topo- 
logically non-trivial configurations of gauge fields known as G-instantons) induce different 
gauge groups on the orbifold fixed planes. For example, it was shown in Refs. ^ and [33 
that one can obtain grand unified gauge groups such as SO (10), SU(5) and the standard 
model gauge group SU(3)c X SU{2) L x U{1)y on the visible brane by appropriate choice of 
G-instanton. Similarly, one can obtain smaller gauge groups on the hidden brane, such as E<j 



and Eq. However, as demonstrated in Refs. 31-33, requiring a physically interesting gauge 



group such as SU(3)c x SU{2) L x U(l)y, along with the requirement that there be three 
families of quarks and leptons, typically leads to the constraint that there must be a certain 
number of M5-branes in the bulk space in order to make the theory anomaly-free. These 
M5-branes are wrapped on holomorphic curves in the Calabi-Yau manifold, and appear as 
3-branes in the five-dimensional effective theory. The five-dimensional effective action for 
non-standard embeddings and bulk space M5-branes was derived in Ref. [31]. The key conclu- 
sion from this body of work is that heterotic M-theory can incorporate the particle content 
and symmetries required for a realistic low-energy effective theory of particle phenomenol- 
ogy. There has been a considerable amount of literature studying both the four- dimensional 
limit of Hof ava-Witten theoryB and heterotic M-theory.il 

A feature of M5-branes which we will utilize is that they are allowed to move along the 
orbifold direction. It is important to note that the M5-brane motion through the bulk is 
not "free" motion. Rather, non-perturbative effects, such as the exchange of virtual open 
supermembranes stretched between a boundary brane and the bulk M5-brane, can produce 
a force between them. The corresponding potential energy can, in principle, be computed 
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from M-theory. Explicit calculations of the supermembrane-induced superpotentials in the 



effective four-dimensional theory have been carried out in Refs. [36] and |3y]. Combined with 
the M5-brane contribution to the Kahler potential presented in Ref. ^8], one can obtain an 
expression for the supermembrane contribution to the M5-brane potential energy. These 
non-perturbative effects cause the bulk brane to move along the extra dimension. In partic- 
ular, it may come into contact with one of the orbifold fixed planes. In this case, the branes 
undergo a "small instanton" phase transition which effectively dissolves the M5-brane, ab- 
sorbing its "data" into the G-instanton.il Such a phase transition at the visible brane may 
change the number of families of quarks and leptons as well as the gauge group. For instance, 
the observable 3-brane may go from having no light families of quarks and leptons prior to 
collision to having three families after the phase transition. 

It is in this five-dimensional world of heterotic M-theory, bounded at the ends of the 
fifth dimension by our visible world and a hidden world, and supporting moving five-branes 
subject to catastrophic family and gauge changing collisions with our visible 3-brane, that 
we propose to find a new theory of the very early universe. 



III. INITIAL CONDITIONS 

All cosmological models, including the hot big bang and the inflationary scenario, rest 
on assumptions about the initial conditions. Despite attempts, no rigorous theory of initial 
conditions yet exists. 

Inflationary theory conventionally assumes that the universe emerges in a high energy 
state of no particular symmetry that is rapidly expanding. If traced backward in time, such 
states possess an initial singularity. This is only one of several fundamental obstacles to con- 
structing a well-defined theory of 'generic' inflationary initial conditions. For a theory based 
on such general and uncertain initial conditions, it is essential that there be a dynamical at- 
tractor mechanism that makes the universe more homogeneous as expansion proceeds, since 
such a mechanism provides hope that the uncertainty associated with the initial conditions 
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is, in the end, irrelevant. Superluminal expansion provides that mechanism. 

The ekpyrotic model is instead built on the assumption that the initial state is quasi- 
static, nearly vacuous and long-lived, with properties dictated by symmetry. So, by con- 
struction, the initial state is special, both physically and mathematically. In this case, while 
a dynamical attractor mechanism may be possible (see below), it is not essential. One can 
envisage the possibility that the initial conditions are simply the result of a selection rule 
dictating maximal symmetry and nearly zero energy. 

Within the context of superstring theory and M-theory, a natural choice with the above 
properties is the BPS (Bogolmon'yi-Prasad-Sommerfeld) state.! The BPS property is already 
required from particle physics in order to have a low-energy, four- dimensional effective action 
with M = 1 supersymmetry, necessary for a realistic phenomenology. For our purposes, the 
BPS state is ideal because, not only is it homogeneous, as one might suppose, but it is 
also flat. That is, the BPS condition links curvature and homogeneity. It requires the two 
boundary branes to be parallel. 

The BPS condition also requires the bulk brane to be nearly stationary. If the bulk 
brane has a small initial velocity, it is free to move along the fifth dimension. Assuming the 
bulk brane to be much lighter than the bounding branes enables us to treat its backreaction 
as a small perturbation on the geometry. Non-perturbative effects can modify this picture 
by introducing a potential for the bulk brane. For example, a bounding brane and bulk 
brane can interact by the exchange of M2-branes wrapped on holomorphic curves, resulting 
in a potential drawing the bulk brane towards our visible brane. In Section V, we shall see 
that the non-perturbative potential plays an important role in determining the spectrum of 
energy density fluctuations following the collision of the bulk brane with our visible brane. 

We do not rule out the possibility of a dynamical attractor mechanism that drives the 
universe towards the BPS state beginning from some more general initial condition. Such 
parallelism would be a natural consequence of all the branes emerging from one parent 
brane. Another appealing possibility is to begin with a configuration consisting of only the 
two bounding 3-branes. The configuration may have some curvature and ripples, but these 
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can be dissipated by radiating excitations tangential to the branes and having them travel 
off to infinity. Then, at some instant, the hidden brane may under go a "small instanton" 
transition which causes a bulk brane to peel off. While little is known about the dynamics 
of this peeling off, it is reasonable to imagine a process similar to bubble nucleation in first 
order phase transitions. We suppose that there is a long-range, attractive, non-perturbative 
potential that draws the bulk brane towards our visible brane, as shown in Figure 2. Very 
close to the hidden brane, there may be a short-range attractive force between the bulk 
brane and the hidden brane due to small-instanton physics (not shown in the Figure). The 
situation is similar to false vacuum decay where the position of the bulk brane along the 
fifth dimension, Y, plays the role of the order parameter or scalar field. Classically, a 
bulk brane attached to the hidden brane is kept there by the energy barrier. However, 
quantum mechanically, it is possible to nucleate a patch of brane for which Y lies on the 
other side of the energy barrier. At the edges of the patch, the brane stretches back towards 
and joins onto the orbifold plane. The nucleated patch would correspond to the minimum 
action tunneling configuration. We conjecture that, as in the case of bubble nucleation, the 
configuration would be one with maximal symmetry. In this case, this configuration would 
correspond to a patch at fixed Y with a spherical boundary along the transverse dimensions. 
The nucleation would appear as the spontaneous appearance of a brane that forms a flat, 
spherical terrace at fixed Y parallel to the bounding branes (i.e., flat). Once nucleated, 
the boundaries of the terrace would spread outwards at the speed of light, analogous to the 
outward expansion of a bubble wall in false vacuum decay. At the same time, the brane 
(the terrace) would travel towards the visible brane due to the non-perturbative potential. 
As we shall see in Section [TV], this motion is very slow (logarithmic with time) so that 
the nucleation rate is essentially instantaneous compared with the time scale of transverse 
motion. In other words, for our purposes, the nucleation process corresponds to a nearly 
infinite brane peeling off almost instantaneously. 

Beginning in an empty, quasi-static state addresses the horizon and flatness problems, but 
one should not underestimate the remaining challenges: how to generate a hot universe, and 
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how to generate perturbations required for large-scale structure. The remarkable feature of 
the ekpyrotic picture, as shown in the forthcoming sections, is that brane collision naturally 
serves both roles. 



V(Y) 




"visible" 
brane 



"hidden" 
brane 



(y=0) (y=R) 

FIG. 2. Sketch of the exponential potential V(Y) = —ve~ maY (the line of zero potential energy 
corresponds to the dotted line). The potential attracts the bulk brane towards the visible brane. 
The force is strongest near the visible brane and tends to zero at large distances. 



In either setup, we begin with a flat bulk brane, either a finite patch or an infinite plane, 
which starts nearly at rest, and a non-perturbative force drawing it towards the visible brane. 
These initial conditions are sufficient to enable our scenario. 



IV. PROPAGATION OF THE BULK BRANE 

A. Moduli Space Actions for Brane- World Gravity 

In this section, we discuss the moduli space approximation which we shall employ 
throughout our analysis. This approximation may be used when there is a continuous 
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family of static solutions of the field equations, of degenerate action. It is the basis for 
much of what is known of the classical and quantum properties of solitons such as magnetic 
monopoles and vortices.0 It is also a powerful tool for cosmology as it neatly approximates 
the five dimensional theory in the regime where the rate of change of the geometry, as mea- 
sured, for example, by the four dimensional Hubble constant H is smaller than the typical 
spatial curvature scale in the static solutions. The moduli are the parameters specifying the 
family of static solutions, 'flat directions' in configuration space along which slow dynamical 
evolution is possible. During such evolution the excitation in other directions is consistently 
small provided those directions are stable and characterized by large oscillatory frequencies. 

The action on moduli space is obtained by substituting the static solutions into the full 
action with the modular parameters represented as space-time dependent moduli fields, Q 1 , 
where / runs over all the moduli fields. If we consider the time dependence first, as we shall 
do for the homogeneous background solutions, the moduli space action takes the form 



where Gjj(Q) is a matrix- valued function of the moduli fields. This is the action for a non- 
relativistic particle moving in a background metric Gu(Q), the metric on moduli space. 
For truly degenerate static solutions, the potential term must be constant, and therefore 
irrelevant to the dynamics. Even if a weak potential V is additionally present, as it shall 
be in our discussion below, the moduli approximation is still valid as long as the dynamical 
evolution consists in the first approximation of an adiabatic progression through the space 
of static solutions. 

In the next section we compute the moduli space action for heterotic M-theory. First, 
however, it is instructive to consider a simpler model which demonstrates similar physical 
effects and is of some interest in its own right. This model, discussed by Randall and Sun- 
drum,0 consists of a five dimensional bulk described by Einstein gravity with a negative 
cosmological constant A, bounded by a pair of branes with tension ±a. The brane tension 
must be fine tuned to the value a = (3|A|/(47rG 5 )) 1//2 in order for static solutions to ex- 




(1) 



17 



ist. Because there are few moduli, the model is simple to analyze and illustrative of some 
important effects. However, as we emphasize below, there are reasons for taking heterotic 
M-theory more seriously as a candidate fundamental description. 

In the Randall- Sundrum model, the static field equations allow a two-parameter family 
of solutions, with metric 

ds 2 = -n 2 dr 2 + a 2 dx 2 + dy 2 , a = e y/L , n = Na, y 1 < y < y , (2) 

and the positive (negative) tension branes located at yo (yi) respectively. L is the Anti-de 
Sitter (AdS) radius given by L 2 = 3/(47t(j5|A|), and N is an arbitrary constant. 

The above solutions are specified by the three moduli N, y and y%. We now allow them 
to be time dependent. The lapse function N(t) is associated with time reparametrization 
invariance and, hence, not a physical degree of freedom. The other two moduli represent the 
proper distance between the branes, but also the time- dependent cosmological 'scale factors' 
ao = e y °l L and a\ = e Vl ^ L on each brane. We shall see in a moment how these combine to give 
four- dimensional gravity with a massless scalar, related to the proper separation between 
the branes. 

To compute the moduli space action it is convenient to change from the coordinates 
in Eq. (El) to coordinates in which the branes are fixed. This is accomplished by setting 
V — (yo ~ Vi)y + Uii an d the branes are now located at y = and 1 respectively. We 
substitute the ansatz Eq. @ into the five dimensional action and integrate over y. All 
potential terms cancel. Since the branes do not move in the y coordinates, the kinetic terms 
arise only from the five dimensional Ricci scalar, which yields the result 

S = [ drd^xN- 1 Q(d\ - a 2 ). (3) 

167TG5 J 

Thus the metric on moduli space is just the 1+1 Minkowski metric, and we infer that moduli 
space in this theory is completely flat. 
We may now change coordinates to 

a = a cosh / a x = a sinh / (4) 
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where Ob ^0 ^1? 

and the proper separation between the branes is just Lln(ao/ai) = 
Lln(coth/). The action (Q) then becomes 

S = -J— [ drd 3 xN^ 6 (-a 2 + a 2 f 2 ) , (5) 

167TG4 J v ' 

where G4 = G^L~ l , just the action for four dimensional gravity coupled to a massless scalar 
field /. This is the 'radion' field.0 Note that the '4d Einstein frame scale factor' a is not 
the scale factor that is seen by matter localized on the branes: such matter sees ao or a%. 
As we shall see, it is perfectly possible for both of the latter scale factors to expand while a 
contracts. 

The general solution to the moduli space theory is easily obtained from ([3]): the scale 
factors evolve linearly in conformal time r, with di = ±do- From the point of view of 
each brane, the motion of the other brane acts as a density of radiation (i.e., allowing 
ai,o — constant). Although the moduli space theory has identical local equations to four- 
dimensional gravity coupled to a massless scalar, the geometrical interpretation is very 
different, so that what is singular from one point of view may be non-singular from the 
other .0 

To obtain an action that might describe the present, hot big bang phase with fixed 
gravitational constant, one might add a potential V(f) which fixes the interbrane separation 



/, causing it to no longer be a free modulus (e.g., see Ref. |42|) . Likewise adding extra fields 
on either brane, one sees that a standard Friedmann constraint is obtained from the variation 
with respect to N. Generally the matter couplings will involve the field /. However, if V(f) 
rises steeply away from its minimum and if its value at the minimum is zero (so that there is 
no vacuum energy contribution), / will not evolve appreciably and four dimensional gravity 
will be accurately reproduced. 

We imagine that, prior to the big bang phase, there is an additional, positive tension 
bulk brane. For static solutions to occur, we require that the three brane tensions sum to 
zero. Positivity of the bulk brane tension imposes that the cosmological constant to its right, 
Ao, be smaller in magnitude than that to its left, Ai, so that the corresponding AdS radii 
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obey L > L\. 

For the three brane case we obtain the moduli space action 

S = y^q~ 5 j dr^xN-'Q (-L a 2 + (L - Li)d| + L^) , (6) 

where Ob is the scale factor on the bulk brane. Again this is remarkably simple, just 
Minkowski space of 2+1 dimensions. Likewise for N parallel branes, all with positive tension 
except the negative tension boundary brane, the system possesses a metric which is that 
for iV-dimensional Minkowski space. Note that the 'masses' appearing in the kinetic terms 
for the boundary branes are the opposite of what one might naively expect. Namely, the 
positive tension boundary brane (hidden brane) has the negative 'mass' — L , whereas the 
negative tension boundary brane (visible brane) has a positive 'mass' L\. The magnitudes 
of these terms are also surprising: the visible brane has the greater magnitude tension, but 
the smaller magnitude moduli space mass. 

Just as for the brane- ant ibrane system, when more branes are present one can change 
variables to those in which the theory resembles four dimensional Einstein gravity coupled 
to massless fields. For three branes, the required change of variables is 

do = a cosh/, 

as = (L /(L — Li))^ asinh/ cos6>, 

ai = (L /Li)2 asinh/sin#, (7) 

and the scalar field kinetic term takes the form a 2 (f 2 + ^ 2 sinh 2 /). The scalar fields live 
on the hyperbolic plane H 2 , and there is a nontrivial Kahler potential. When we introduce 
potentials, we must do so in a manner which respects four dimensional general coordinate 
invariance. This restricts the form to 

AS = - J dTd 3 xNa 4 V{f, 6). (8) 

We shall assume that the system starts out nearly static, and that the potential energy 
is always negative, so that we are clearly not using inflation to drive expansion. We further 
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assume that the interaction potential draws the bulk brane away from the hidden brane 
towards the visible brane. The original variables in Eq. (|6]) provide some insight into what 
happens. We consider an interaction between ao and a# causing the latter to decrease. 
But since ao has a negative 'mass' it is actually pushed in the same direction and thus 
contracts. Similarly an interaction potential between ao and a\ can have the opposite effect: 
causing both ao and a\ to expand. Figure (|f) shows an example, with the bulk brane being 
pushed across the gap, and the visible brane being attracted towards it. The hidden brane 
actually 'bounces' (this is barely visible in the Figure) due to a competition between the two 
effects, so that by the collision between bulk and visible branes, both outer boundary branes 
are actually expanding; that is, a and ai are both increasing. At first sight this appears 
inconsistent with the four dimensional point of view: if the system starts out static, and with 
all potentials negative, then the 4d Einstein frame scale factor a must contract throughout. 
The two points of view are consistent because di contributes negatively to d. Since a\ is 
expanding rapidly compared to a , a is indeed contracting, as shown in the Figure. Of 
course what is happening physically is that the fifth dimension is collapsing, a well known 
hazard of Kaluza-Klein cosmology. Here, the inter-brane separation is decreasing while the 
scale factors seen by matter on the branes are expanding. 
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FIG. 3. Evolution of scale factors in a three-brane system where a bulk brane is drawn across 
from the hidden brane to the visible brane. The solid lines show (from top to bottom) the scale 
factors on the hidden brane (ao), the bulk brane (as) and the visible brane (ai). At collision 
between the bulk brane and the visible brane, both boundary branes are expanding. The evolution 
of the 4d Einstein frame scale factor a is also shown as a dashed line: it contracts throughout. 

What happens when the bulk brane meets the boundary? A matching condition is needed 
to determine the resulting cosmology. The initial state is specified by the scale factors ao, 
as and a\ and their time derivatives, the final state by ao, a\ and their time derivatives, 
plus the coordinates and momenta of any excitations produced in the collision. 

One expects that the scale factors a, a\ and ao should be continuous: likewise do would 
be expected to be continuous if no bulk brane hits the hidden brane, and if no constraint 
on the size of the extra dimension is imposed. However, di cannot be continuous since the 
bulk brane imparts some momentum on the visible brane. The momentum conservation 
condition can be expressed as 
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where i and / subscripts label initial and final quantities. Conservation of momentum applies 
if the forces derive from a potential V which is short-ranged and translation invariant and if 
there is no other entity that carries momentum after collision. The above conditions imply 
the continuity of the '4d Einstein frame scale factor' a and its time derivative. 

This matching condition, if correct, would pose a serious problem for our scenario since 
it implies that the four dimensional scale factor a is contracting after collision, at the be- 
ginning of the hot big bang phase. Specifically, the matching condition suggests a simple 
continuation of the motion in Figure ||| after collision in which both branes are expanding but 
a is decreasing because the two branes are approaching one another and the fifth dimension 
is collapsing. 

In the M-theory models which we consider in this paper, we shall simply impose a 
constraint on moduli space which ensures that the distance between the visible and hidden 
branes becomes fixed after collision, as required to converge to the Hofava-Witten picture. 
The constraint corresponds to fixing a^/ai after collision, forcing a discontinuity in both 
di and do- In this case the momentum matching condition (continuity of d) yields rather 
paradoxical behavior in which both di and do reverse after collision, so the universe collapses. 
This does not seem physically plausible, especially when matter is produced at collision, and 
the expansion of the matter would then have to be reversed as the size of the fifth dimension 
became fixed. More plausible is that d is also discontinuous: the branes collide, their 
separation becomes fixed, and the pair continue in the same direction of motion (expansion) 
as before collision. Here we simply wish to flag this issue as one that we have not resolved 
in the M-theory models considered here: more work is needed to do so. 

We have identified at least one mechanism for avoiding contraction or collision while 
still remaining within the moduli space approximation. We have constructed models for 
branes in AdS employing 'non-minimal' corrections to the kinetic terms of a^o, which are 
allowed by four dimensional general coordinate invariance. These non-minimal kinetic terms 
both stabilize the size of the extra dimension and allow final expansion from static initial 
conditions, with negative potentials.^ 
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Another possibility is that the moduli space approximation break-down at collision (it 
must break down, since radiation is produced) leads to the release of radiation into the 
bulk. This is prohibited by planar symmetry in the AdS example, but is possible in the 
more general M-theory context. Radiation emitted into the bulk contributes to the pressure 
T 55 , which, from the G 55 Einstein equations, acts to decelerate a±. The emitted radiation is 
redshifted as it crosses the bulk, so is likely to have less effect on the hidden brane if it is 
absorbed there. The net result would be a slowing of a±, causing the effective scale factor a 
to increase. 

We do not want to understate the challenge of obtaining a final expanding universe with 
stabilized fifth dimension. In a conventional four-dimensional theory (Einstein gravity plus 
scalar fields) it would simply be impossible to start from zero energy and, through evolution 
involving negative potentials, obtain a final expanding universe. Our point is that brane 
world scenarios offer ways around this 'no-go theorem', which we have just begun to explore. 

The AdS examples we have discussed are instructive in that they are easier to ana- 
lyze than the full M-theory case. However, as mentioned above, it is unlikely that these 
model theories are quantum mechanically consistent. The most obvious problem is that fine 
tuning is needed to balance the brane tension against the cosmological term. Without this 
balance, no static solutions are possible. Computing the quantum corrections may in fact be 
impossible since, in the thin-brane limit, these are generally infinite and non-renormalizable. 

Therefore, for the remainder of the paper, we turn to analogous examples in heterotic 
M-theory, which is more complex but has other advantages. The branes in this theory are 
BPS states, protected from quantum corrections by supersymmetry. Their tensions are fixed 
by exact quantum mechanical symmetries and there is no fine tuning problem analogous to 
that present in the Randall-Sundrum models. 
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B. The Background BPS Solution in Heterotic M-theory 



The five- dimensional effective action of heterotic M-theory was derived in Refs. y and 30. 
Its field content includes a myriad of moduli, most of which will be assumed frozen in this 
paper. We shall, therefore, use a simplified action describing gravity g^s, the universal 
"breathing" modulus of the Calabi-Yau three-fold (ft, a four-form gauge field A^c with field 
strength T = dA and a single bulk M5-brane. It is given by 



s = ^f d^^-g[n-Ud<pf 

L JMf. \ 4 



I m 5 V 2 2 5!/ 

3 „ / fiunX \ 

-3 £ a t Mt j Mf d% l} ^f%e^ - —A^d.X^X^X^xU , (10) 

where 7, 5, e, £ = 0, . . . , 4, /1, z/, . . . = 0, . . . , 3. The space-time is a five-dimensional manifold 
A4§ with coordinates x 1 . The four-dimensional manifolds .M4 , i = 1,2,3 are the visible, 
hidden, and bulk branes respectively, and have internal coordinates £^ and tension OjMf . 
Note that cVj has dimension of mass. If we denote a% = —a, = a — j3, and 0:3 = (3, 
then the visible brane has tension — aMf , the hidden brane (a — /3)M|, and the bulk brane 
/3M|. It is straightforward to show that the tension of the bulk brane, /3M|, must always be 
positive. Furthermore, one can easily deduce that the tension on the visible brane, — oMf , 
can be either positive or negative. The ekpyrotic scenario can be applied, in principle, to 
any such vacua. In this paper, for specificity, we will always take a > 0, so that the tension 
on the visible brane is negative. Furthermore, we will choose f3 such that a — (3 > 0, that is, 
the tension of the hidden brane is positive. The tensor hfy is the induced metric (and 

'(i)V«i(i)> 

111 uiiiu wuius, Xfctefl, 

into M5. 

The BPS solution of Lukas, Ovrut, and Waldrami is then given by[] 



its determinant) on . The functions XZJ^jL) are the coordinates in Ai 5 of a point on 
AI4 with coordinates In other words, XZJ^) describe the embedding of the branes 



1 We have changed the notation used in Ref. |6| by replacing their H(y) with D(y). In this paper, 
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ds 2 = D(y)(-N 2 dr 2 + A 2 dx 2 ) + B 2 D\y)dy 2 
= BD\y) 

-foi23Y = -aA^B-'D- 2 ^) for y < Y 

= -(a - (3)A 3 NB~ 1 D~ 2 (y) for y > Y, (11) 

where 

D(y) = ay + C for y < Y 

= (a- fl)y + C + I3Y for y > Y, (12) 

and A,B,C,N and Y are constants. Note that A,B,C,N are dimensionless and Y has the 
dimension of length. The visible and hidden boundary branes are located at y = and 
y = R, respectively, and the bulk brane is located at y = Y, < Y < R. We assume that 
C > so that the curvature singularity at D = does not fall between the boundary branes. 
Note that y = lies in the region of smaller volume while y = R lies in the region of larger 
volume. 

Finally, note that inserting the solution of the four-form equation of motion into Eq. ( [X0|) 
yields precisely the bulk action given in Ref. ^ with charge —a in the interval < y < Y 
and charge — a + (3 in the interval Y < y < R. The formulation of the action Eq. ( |T0| ) using 
the four-form A is particularly useful when the theory contains bulk branes, as is the case 
in ekpyrotic theory. 



C. The Moduli Space Action of Heterotic M-theory 

As in Section [IV A|, we shall use the moduli space approximation to study the 



dynamics of heterotic M-theory with a bulk brane. The static BPS solution in- 
volves five constants A, B, C, N, and Y. These now become the moduli fields, Q 1 = 



we will use the symbol H to denote the Hubble parameter. Furthermore, comparing with their 
notation, we have rescaled a by a factor of \/2/3 and have defined = V . 
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(A(x, r), B(x, r), C(x,t), N(x,t), Y(x,t)). In the limit of homogeneity and isotropy, the 



moduli fields are functions of time only. Substituting the static ansatz ( |TT]) into the ac- 
tion (|I0D, and integrating over y, we obtain the moduli space action S mo d with Lagrangian 
density 



£ 



mod 



G IJ (Q)Q I Q J -V(Q) = C bulk + £ 



Pi 



(13) 
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and 
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[(aY + C) m+1 - C m+1 ] 

[((a-(3)R + C + (3Y) m+1 



{aY + C) 



m+l] 



(15) 



Note that I m has the dimension of length. We see from Eq. (|H"D that the Lagrangian of the 
4d effective theory is the sum of two parts, £buik and Cp. The first contribution, £buik, comes 
from the bulk part of the five-dimensional action, whereas the second contribution, £@, is 
the Lagrangian of the bulk brane. Note that we have added by hand a potential V(Y) in Cp 
which is meant to describe non-perturbative interactions between the bulk brane and the 
boundary branes.&il The actions of the two boundary branes, which are at fixed values of 
y, do not contribute to the 4d effective action. Their contribution is canceled by bulk terms 
upon integration over y. This cancellation is crucial, since it yields a 4d effective theory 
with no potentials for A, B, C, or N, thereby confirming that these fields are truly moduli 
of the theory. 
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Eq. fll4|) is analogous to the action for gravity with scale factor A coupled to scalar 
fields. Since the overall factor of BI 3 M 3 will generically be time- dependent before collision, 
it follows that the scalar fields B, C, and Y are non- minimally coupled to gravity. In order 
to match onto a theory with fixed Newton's constant G4, we shall impose the condition that 
this factor become constant after collision. Hence, after collision, we can identify the 4d 
effective Planck mass as 



where ^3 is evaluated at Y = and for B,C at the moment of collision. Note that in the 
limit (a - C C, this expression agrees with the 4d Planck mass identified in Ref. [15. 

At this point, one can define a new scale factor a = A{BI^M^) 1 / 2 , as well as n = 
Ni^BI^M^) 1 / 2 . This has the effect of removing off-diagonal terms in the moduli space metric 
that couple the field A to the other variables. In these new variables, the bulk Lagrangian 
becomes 



3a 3 M 2 fa\ 2 1 f BY IB I 



^{(8-£)Mw-S)^8-£K 



In this form, it is clear that a is an scale factor analogous to the variable a of the previous 
section. The moduli B, C and Y behave effectively as scalar fields, albeit with nontrivial 
kinetic terms. 

In analogy with the example in Section IV A and the discussion following Eq. (||), we 
now impose two constraints consistent with four dimensional covariance. These reduce the 
moduli degrees of freedom and simplify the system. Namely, we shall impose that 

B = constant 

C = constant. (18) 

At the moment of collision, the modulus Y disappears from the theory. The above conditions 
then imply that the distance between the boundary branes as well as the volume of the 
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Calabi-Yau three-fold become fixed. This is necessary if we want to match onto a theory with 
fixed gravitational and gauge coupling constants. For instance, since ^3 becomes constant 
at collision, it follows from Eq. ([16]) that M p i freezes at that point. 

If we impose Eq. (|IBD , and if we further assume that the tension of the bulk brane is 
small compared to that of the boundary branes, that is, (3 <C a (which allows us to neglect 
the correction to the kinetic term for Y coming from Cbuik) , then the full Lagrangian reduces 
to 

£ _ 3M 5 V f /d\ 2 | g 
n \ \aj I 3 

This action describes a scalar field Y minimally coupled to a gravitational background with 
scale factor a. Note that the gravitational coupling constant associated with a is simply M 5 . 

A repeating theme of this paper is that the effective scale factor (and the associated 
Hubble parameter) can be defined several different ways involving different combinations of 
moduli fields, depending on the physical question being addressed. For example, looking 
ahead, we will show that a is not the same as the effective scale factor for an observer on a 
hypersurface of constant y or the effective scale factor relevant for describing fluctuations in 
Y. 



l D{ YfY 2 -n 2 ^Q- 
2 1 ; BhM, 



(19) 



D. Equations of Motion, the Ekpyrotic Temperature, and the Horizon Problem 



From the moduli space action obtained in Section |IV Q , we can find the equation of 
motion for the bulk brane, described by Y, in the limit that the bulk brane tension is small, 
/? < a. We will use this equation to compute the ekpyrotic temperature, the temperature 
immediately after the branes collide and the universe bursts into the big bang phase. We 
will then compare the Hubble radius at the beginning of the big bang phase to the causal 
horizon distance estimated by computing the time it takes for the bulk brane to traverse the 
fifth dimension. 

Recall that for the static BPS solution, A and N are constants. Without loss of generality, 
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we can choose A = N = 1 in the BPS limit. Variation of Eq. ([19]) with respect to n yields 
the Friedmann equation which, setting n = a, is given by 

H^i^X = J^i\D(YYY* + V{Y)). (20) 



Here we have introduced H to denote the Hubble constant associated with a. Since the 
gravitational coupling constant associated with a is M 5 , we can identify from Eq. (|20f) the 
energy density of the bulk brane 

P( s = 3M*H 2 . (21) 
The equation of motion for a yields the second equation of FRW cosmology 



^ (^-D(y) 2 y 2 - 2V(y)J • (22) 

Finally, we can express the equation of motion for Y in a simple way by defining $ such 
that * = (£>(y)/(/ 3 M 5 ) 1 / 2 )y. Once again making the gauge choice n = a, one finds that 
\& satisfies 

AQ a -^ + ^(*))^-3(A)^, (23 ) 

where V e /j = V(Y)/(B(I 3 M 5 ) 2 ). In this form, Eq. (f23|) looks like the equation of motion 
for a scalar field in a cosmological background. It is, therefore, simple to analyze. The left 
hand side is the time derivative of the total energy density associated with the motion of Y. 
The right hand side either decreases or increases due to the cosmic evolution described by 
the scale factor a. Hence, as regards the kinetic energy of the bulk brane, a is the relevant 
scale factor. 

Perturbing around the BPS limit, we have A = l + 0((3/a) and -Z3 = 7^ +(D((3/a) 7 where 
is the value of ^3 when (3 — and is time-independent. Therefore, a = M 5 ) 1//2 + 

0(P/a). If V(Y) < 0, then all contributions on the right hand side of Eq. ( |22| ) cause a to 
be negative. Hence, beginning from a static initial condition (d = 0), a will contract. If a 
contracts, then, from Eq. (^), the total energy of ^ grows with time. Thus, as long as V 
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is negative semi-definite, the field ^ gains energy as the bulk brane proceeds through the 
extra dimension. 

At the moment of collision, the modulus Y effectively goes away and new moduli describ- 
ing the small instanton transition and new vector bundle take its place. In particular, the 
potential for Y matches onto the potential for these new moduli. While V(Y) will generi- 
cally be negative up to the moment of collision, we shall assume that, once the other moduli 
are excited during the small instanton transition, the potential rises back up to zero in the 
internal space of the new moduli such as to leave the universe with vanishing cosmological 
constant. We base this assumption on the notion that both the initial and final states consist 
of two branes in a BPS vacuum with V = 0. That is, if we were to adiabatically detach the 
bulk brane from the hidden brane and, then, transport and attach it to the visible brane, 
both the initial and final states would consist of two branes only in a BPS configuration 
with cosmological constant zero. 

To mimic the effect of the other moduli during the small instanton transition, we shall 
assume that V(Y) is negative and approaches zero as Y — > 0. Note that if a were constant, 
the kinetic energy and the total energy, p@, would go to zero as V — > and as the branes 
collide. However, from Eq. Q23D, we see that \& (and therefore Y) has extra kinetic energy 
as V — > due to the gravitational blue shift effect caused by a contracting scale factor 
a. We assume that the extra kinetic energy (equal to p@ at collision) is converted into 
excitations of light degrees of freedom, at which point the radiation-dominated era begins. 
The temperature after the conflagration that arises from the brane collision is referred to as 
the "ekpyrotic temperature," analogous to the reheat temperature after inflation. 

Conceivably, some fraction of the extra kinetic energy is converted into thermal excita- 
tions, some into coherent motion of moduli fields, and some, perhaps, into bulk excitations 
(gravitons). If coherent motion is associated with massless degrees of freedom, the associ- 
ated kinetic energy redshifts away faster than radiation. If associated degrees of freedom 
are massive, they can ultimately decay into radiation. Neither case is problematic. For 
simplicity, we shall just assume that the kinetic energy at collision is converted entirely into 
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radiation with an efficiency of order unity. 

The collision energy can be computed by integrating the equation of motion, Eq. (f23|), 
which is expressed explicitly in terms of the time derivative of the total bulk brane energy. 
However, it is somewhat simpler to first compute the Hubble parameter H c upon collision 
using Eq. Q2"2"|), and then to substitute in Eq. fl2"T| ) to obtain the collision energy pp. (Note 
that the subscript c denotes that the quantity is evaluated at the moment of collision.) 
Eq. ( f22|) can be rewritten as 

- « -J- {\d{Y?Y 2 + V(Y) - WOO) « -J- (-3V(y)), (24) 

where we have used the fact that, to leading order in P/a, 

-D 2 {Y)Y 2 + V{Y) = E « 0. (25) 

Note that E is the total energy of the bulk brane which is assumed small compared to the 
energy gained from gravity. We can then integrate Eq. (pll) to obtain 

n r Y=0 R A T r Y=R 

- = -3 / f {~V{Y))%dY = -J? / D(-V)%, (26) 

where we have made use of the fact that I3 = J^ ^ = constant to leading order in (3/ a. Since 
(BI^Ms) 1 / 2 , we obtain 



a 2 



5 



C R D(-V) l/2 dY. (27) 
jy=o 



V2B(I 3 M 5 ) 3 / 2 Jy=o 

Now that we have an expression for the Hubble parameter H c , we can substitute in Eq. fl2~ID 
and find 



The corresponding ekpyrotic temperature is then 



/ D(-y) 1 / 2 dy 

JY=0 



1/2 



(29) 



M pl \MI) (2 J B) 1 /4(j 3 M 5 )3/4 
(N.B. The identification of energy density or effective Planck mass may vary under Weyl 
transformation, but the ratio T/M p i is invariant.) For instance, consider a potential of the 
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form V(Y) 



= —ve 



, where v and m are positive, dimensionless constants. Since non 



perturbative potentials derived from string and M-theory are generically of exponential form 
(for motivation, see the discussion under Eq. (|64|)), this potential will be a standard example 
throughout. In that case, the temperature is calculated, using Eq. (|29|), to be 



where we have used Eq. flT6l). As an example, we might suppose a = 2000M 5 , (3 = M 5 , 
B = KT 4 , C = 1000, R = Mg" 1 , v ~ 10~ 8 , and m = 0.1, all plausible values. This gives 
M5 = 10~ 2 M P ; and produces an ekpyrotic temperature of 10 11 GeV. Note that, with these 
parameters, the magnitude of the potential energy density for Y is (10 _6 M p /) 4 at collision. 
Thus, the typical energy scale for the potential is 10 13 GeV. Later, we will see that these same 
parameters produce an acceptable fluctuation amplitude. We want to emphasize, however, 
that there is a very wide range of parameters that lead to acceptable cosmological scenarios. 

Let us now turn our attention to the homogeneity problem. We have argued that the 
universe begins in a BPS state, which is homogeneous. This condition is stronger than 
needed to solve the homogeneity problem. It suffices that the universe be homogeneous on 
scales smaller than the causal, particle horizon. Let (— r) to t denote the time taken by the 
bulk brane to travel from the hidden to the visible brane. By integrating Eq. (p5|), we find 
that the comoving time is 



The horizon distance dnoRi as measured by an observer on the visible brane, is the elapsed 
comoving time at collision (Y = 0) times the scale factor, D 1 ^ 2 (y = 0) = C 1 / 2 . We find, for 
an exponential potential of the above form, that 




(30) 




(31) 




(32) 



On the other hand, the Hubble radius at collision is obtained from Eq. (|27| ) 




(33) 
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The causal horizon problem is solved if the particle horizon at collision satisfies 
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The condition is easily satisfied for the values of parameters mentioned above, where 
maR/2 ~ 10 2 , the equivalent of 100 e-folds of hyperexpansion in an inflationary model. 

E. Cosmological Evolution for an Observer at Fixed y 

We have seen that the scale factor relevant to describing the equation of motion for Y is 
a, a particular combination of moduli fields. Furthermore, if we assume nearly BPS initial 
conditions (that is, vanishing potential and kinetic energy), then a is a decreasing function 
of time. Hence, the scalar field Y evolves as if the universe is contracting. 

However, as pointed out previously, an important feature is that the effective scale fac- 
tor for other physical quantities depends on other combinations of moduli fields. In this 
subsection, we shall derive the cosmological evolution as seen by an observer living on a 
hypersurface of constant y (for example, an observer on the visible brane). We will find that 
the scale factor seen by such family of observers is different than a. In particular, we find 
that any such observer sees the universe expanding before the bulk brane collides with the 
visible brane. 

Consider, for concreteness, an observer living on the visible brane. The induced metric 
on that hypersurface is obtained from Eq. (|TT|) 

dsl =Q = a C (-dr 2 + dx 2 ) = al(-dr 2 + dx 2 ). (35) 
BI 3 M 5 

The rate of change of the induced scale factor a\ can be written as 



ai d Ju^r d I 2b ^-WiX) 

— = -- 3/3— Y = - + 3(3— , 36 

ax a i 3 a 1 3 1J{Y) 



where we have used Eq. (p5|). In this way, we have expressed di as the sum of two contribu- 
tions: the first contribution is given by a and tends to make a\ contract; the second term, 
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coming from J3, is positive and tends to make ai expand. To determine which of these two 
terms dominates at the moment of collision, we note from Eq. ( pip that 

- ~ {-V(Y)) 1 ' 2 . (37) 

Therefore, both terms on the right hand side of Eq. (^) are proportional to (— ViY)) 1 ^ 2 . 
However, the coefficients are different functions of a and C . For the case of the exponential 
potential, for instance, one finds that reasonable values of a, C, and m (such as those given 
at the end of Section |IVD|) result in the I3 term being larger than the a term in Eq. (j36[). 
The net effect is to make a\ grow with time; that is, an observer at y = sees the universe 



expanding. This is in agreement with the results for the AdS case presented in Section [TV A . 

For other hypersurfaces of constant y, a similar story holds. The induced scale factor 
is the product of a which decreases and a function of Y which increases. Once again, 
reasonable values of the parameters result in expanding hypersurfaces at the time of collision. 
In particular, do, the scale factor on the y = R hidden brane, is expanding at collision, in 
agreement with the AdS results. 

To summarize, for shallow potentials (for example, an exponential potential) and rea- 
sonable values of the parameters, we have seen that both ai and are expanding at the 
moment of collision. On the other hand, a is contracting. Since this agrees qualitatively 
with the AdS case (see Fig. 3), the discussion at the end of Section |IV A| concerning the 
matching condition at collision and the subsequent expansion of the universe applies also to 
heterotic M-theory. 



F. Cosmological Evolution for an Observer on the Bulk Brane 

So far, we have adopted the point of view of the low-energy four-dimensional effective 
action. It is sometimes useful, for instance in calculating the fluctuation spectrum, to adopt 
the point of view of an observer living on the bulk brane. The motion of the bulk brane 
through the curved space-time induces an FRW evolution on its worldvolume. In terms of 
conformal time 77 on the bulk brane, the induced metric is given to lowest order in Y by 
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ds\ = D(Y( v ))(-dr] 2 + dx 2 ). (38) 

From the form of the metric in Eq. (|38|), we see that the scale factor a#(?7) describing the 
FRW evolution on the brane is simply given by clb^i]) = D 1 ' 2 ^). Since the bulk brane 
moves from a region of larger D (location of the hidden brane) to a region of smaller D 
(location of the visible brane), an observer on the bulk brane sees a contracting universe (as 
opposed to an observer on the visible brane who sees expansion). Finally, we note that for 
non-relativistic motion, one has rj « r, where r is global conformal time. 

G. Summary of homogeneous propagation of the bulk brane 

In this section, we have described the spatially homogeneous propagation of the bulk 
brane in terms of the evolution of three different quantities: the scale factor as felt by the 
modulus Y of the bulk brane, the scale factor for an observer living on a hypersurface of 
constant y, and the scale factor for an observer on the bulk brane. For nearly BPS initial 
conditions, we have seen that the scale factor that appears in the equation of motion for 
the bulk brane, namely a, decreases with time. This means that there is a gravitational 
blue shift effect that increases the kinetic energy of the bulk brane. This added energy, we 
propose, is converted to radiation and matter upon collision. Furthermore, we have shown 
that observers at fixed y generically see an expanding universe. Finally, an observer on the 
bulk brane sees a contracting universe. This is a simple geometrical consequence of the fact 
that the bulk brane travels from a region of smaller curvature to a region of larger curvature. 

We should mention that both the modulus Y and observers at fixed y see scale factors 
which are slowly- varying in time, in fact almost constant. This is because both variations 
are due to the back-reaction of the bulk brane onto the geometry, an effect which is of order 
(3 /a. On the other hand, the cosmological evolution felt by an observer on the bulk brane 
is faster, although only by a logarithmic factor. No superluminal expansion is taking place 
from the point-of-view of any observer. Rather, what characterizes the ekpyrotic scenario 
is that all motion and expansion is taking place exceedingly slowly for an exceedingly long 
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period of time. 



V. SPECTRUM OF FLUCTUATIONS 



In this section, we show how the ekpyrotic scenario can produce a nearly scale-invariant 
spectrum of adiabatic, gaussian, scalar (energy density) perturbations that may account for 
the anisotropy of the cosmic microwave background and seed large-scale structure formation. 
The density perturbations are caused by ripples in the bulk brane which are generated by 
quantum fluctuations as the brane traverses the bulk. The ripples result in 3d spatial varia- 
tions in the time of collision and thermalization, and, consequently, they induce temperature 
fluctuations in the hot big bang phase. 

Because both the ekpyrotic scenario and inflationary cosmology rely on quantum fluc- 
tuations to generate adiabatic perturbations, the calculational formalism for predicting the 
perturbation spectrum and many of the equations are remarkably similar. One difference is 
that inflation entails superluminal expansion and the ekpyrotic scenario does not. For the 
ekpyrotic scenario, the fluctuations are generated as the bulk brane moves slowly through 
the bulk. For the examples considered here, the motion is in the direction in which the warp 
factor is shrinking. Because of the shrinking warp factor, the Hubble radius for an observer 
on the brane is decreasing. The effect of a decreasing Hubble radius is to make the spectrum 
blue. In inflation, the Hubble radius is expanding in the 4d space-time, and, consequently, 
the spectrum is typically red. 

Here, we give an abbreviated version of the derivation that emphasizes the similarities and 
differences from the inflationary case. For this purpose, we adapt the "time-delay" approach 
introduced by Guth and Pi for the case of inflation to the colliding brane picture. Eil This 
approach has the advantage that it is relatively simple and intuitive. Experts are aware 
that this approach is inexact and non-rigorous@ and, hence, might question the reliability. 
The more cumbersome and less intuitive gauge invariant approach introduced by Bardeen, 



Steinhardt and 




preferable since it is 
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rigorous and applies to a wider range of models. We have developed the analogue of the 
gauge invariant approach for the colliding branes, and we find that the time-delay approach 
does match for the cases we consider. We will present the gauge invariant formalism for the 
ekpyrotic scenario and a fully detailed analysis in a separate publication.^ 

In the first subsection, we shall assume that the ripples in the bulk brane have already 
been generated and begin our computation just as the bulk brane collides with our visible 
brane. The bulk brane has position Y(t,x) = Yq(t) + SY(r,x), where Yq(t) is the average 
position of the brane along the bulk (y) direction and SY represents the small ripples. Our 
goal is to adapt the time-delay formalism to compute how the ripples translate into density 
fluctuations in the hot big bang phase. In the second subsection, we shall discuss how 8Y is 
set by quantum fluctuations and the general conditions under which the fluctuation spectrum 
will be nearly scale-invariant. Then, we will present specific models that satisfy the scale- 
invariant conditions and discuss general model-building principles and constraints. In the 
final subsection, we will present the computation for the case of tensor (gravitational wave) 
perturbations and show that the spectrum is tilted strongly towards the blue, a prediction 
that differs significantly from inflationary models. 



A. From brane ripples to density fluctuations 

The fluctuations 5Y result in variations in the time of collision (St) that depend on 
the position along the bulk surface, x. In this sense, the bulk brane position Y plays a 
role analogous to the inflaton and the fluctuations SY play a role similar to inflaton 
fluctuations 5<fi. The time-delay formalism applies under the assumption that the time 
delay is independent of time when the perturbations are well outside the horizon; that is, 
St = 8t(x). The formalism, then, allows one to compute how St(x) converts into a density 
perturbation amplitude. 

In the de Sitter limit, one has St ~ S(f>/4>, where the fluctuations 8<j) are time-independent 
and <fi is also time- independent. Hence, the assumption of the time-delay formalism is 
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satisfied. In the ekpyrotic scenario, St ~ SY/Y. Note that SY is time-dependent, and 
so is Y. However, under circumstances to be discussed later in this section, the two have 
the same time-dependence and, consequently, St is time-independent, as required. (The 
time-independent condition is only approximate in both scenarios. The weakness of the 
time-delay approach is that it cannot be simply generalized to the time-dependent case; 
corrections must be computed using a gauge invariant formulation.il) 

We shall assume that the stress-energy tensor after collision is that of an ideal fluid with 
pressure P and energy density p = 3P 

T „u = pg „u + ( P + p) u » u ^ (39) 

where w M is the velocity of the fluid normalized to u 2 = — 1. 

The perturbations can be characterized by the OlsonE! variables, S and div X, defined 

by 

S = -l + 3 P -(/^V^)- 2 

divX = V^Vvp), (40) 

where W ee g^ v + u^u v . The calculation then proceeds in two steps. First, we find the 
value of S and divX at the moment of collision. Secondly, we calculate the time-evolution 
of S in a radiation-dominated universe. 

If the average collision time is r = 0, then r' ee t — St(x) = is the time when collision 
occurs at position x. We have 

divX (r' = 0) = Qa- 2 H c H c ^ST(x). (41) 
It is useful to define a dimensionless time variable x by 

*=U T ' + iid- (42) 
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When x < 1 (x > 1), the mode is outside (inside) the sound horizon. Then, as shown by 
Olson, the density perturbation with wavenumber k, 5k, is related to the Fourier mode S& 
via 

ii^ = s - (43) 

where Sk satisfies the evolution equation 

*1^- 2 ^ + < 2 + * 2 > St = a < 44) 

The solution of the Sk equation can be written as 

Sk = C\x sin x + C2X cos x. (45) 

To fix the coefficients C\ and C2, we use the initial conditions given in Eq. ([|l|). In terms 
of x, the collision time is x c = kj (\^3a c H c ), and the conditions at collision read 

S k = 2x 2 c a c H c Ar(k) 

= 2 1 1 - A) s " (46) 

where Ar(/c) = k 3 / 2 5r(k) / (2n) 3j/2 . Using these conditions, one finds that the coefficients C\ 
and C2 are given by 

The modes of interest lie far outside the horizon at the time of collision, that is, x c <C 1. 
Thus, when x ^> 1 (when the mode comes back inside the horizon), the second term on the 
right hand side of Eq. (|45|) is suppressed by a factor of x c and is therefore negligible. (This 
is the "decaying" mode.) Using this fact in integrating Eq. (f4*3D, one obtains 

|4| = 4a c # c (l - \Ar(k)\. (48) 
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We see that 5k is the product of two factors: the factor a c H c \ Ar(k) | accounts for the fact that 
different regions of space heat up and therefore begin to redshift at different times, while the 
factor in parentheses depends on H c and describes how the change in the Hubble parameter 
during collision affects the fluctuations. For inflationary models near the de Sitter limit, 
H c — > 0, and so 5k is directly related to the time delay 5r(k). For the ekpyrotic model, the 
scale factor a(r) is of the form const. + log r and H/aH 2 is approximately constant. Hence, 
once again, 5k is directly related to the time delay. In both scenarios, Eq. (f|8]) agrees with 
the exact gauge invariant calculation of density perturbations except for small corrections 
to the prefactor. 

For example, consider the exponential potential V = —ve~ maY . In that case, Eq. (f48|) 
yields 



. . . Am 2 ay2v . . , . , tnS 

141 " ~mC + 2~ (49) 

where we have used Eq. (pTD, and where we have assumed that H C J (a c H 2 ) ^> 1, a reasonable 
approximation for the values of a, C, and m of interest. 

B. Prom quantum fluctuations to brane ripples 

Eq. (|48| ) expresses the density perturbation in terms of the time delay at the time of 
collision, Ar(fc). In this section, we compute the spectrum of quantum fluctuations of the 
brane 5Yk and use the result to compute the time delay, Ar(k). 

1. The scalar fluctuation equation 

For the calculation of quantum fluctuations, it is sufficient to work at the lowest order 
in f3/a. Without loss of generality, we can therefore set A — N — 1. In that case, the bulk 
brane Lagrangian is given by 



Cp = 3[3M*B 



l -D(Y) 2 ^d,Yd v Y - V(Y) 
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(50) 



Note that this agrees with Cp given in Eq. (|T4] ) when we set A = N = 1 and spatial gradients 
of Y to zero. Let us first consider the spatially homogeneous motion of the brane which will 
be described by Yq(t). (The subscript "0" emphasizes that we want to think of Yo as the 
background motion.) It is governed by the following equation of motion 

1 



-D(Y ) 2 Y + V(Y ) = E, 



(51) 



where E is a constant. Eq. (^l|) is, of course, simply the statement that the energy E of the 
bulk brane is conserved to this order in (3 /a. Since we have chosen the visible brane to lie 
at y = and the hidden universe to lie at y = R, we focus on the branch Y < in which 
case the bulk brane moves towards the visible brane. The solution to Eq. ([51]) is then given 
by 

D(Y')dY' 



(52) 



'2{E-V{Y')) 

with r < 0, and with the collision occurring at r = 0. 

Let us now consider fluctuations around the background solution Y (t). Namely, if 
Y = Yq(t) + 5Y(t, x), with 5Y(t, x) Yq(t), we can expand the action to quadratic order 
in 5Y 



C 



1 



flue 



D 2 [-5Y 2 + (d(5Y)) 2 ] + 



a 2 D^(V -E)-aD^-^ + 



dY ' 2dY 2 



(SY) 2 . (53) 



where we have used Eq. (|5T|), and where we have introduced D = D(Y ) and Vq = V(Yq) 
for simplicity. 

The key relation is the fluctuation equation as derived from the action (|53|) 



x 



d 2 

2 a Jk 

dx 2 



''pert 



_ Qpert 



X 



x=\k\(-r), 



(54) 



where /r = Dq • SYz and where a pert is defined by 

dV „ d?V n 



dpert 
Qpert 



n~ 3 



oc- 



dYn 



-D, 



dY 2 



(55) 
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The fluctuation equation, Eq. fl54|), can be compared with the corresponding equation 
for the perturbations of a scalar field with no potential and minimally coupled to an FRW 
background with scale factor a(r) 



% + 2~<% + k 2 5^ = 0. (56) 
Defining f^ = a ■ <50£, Eq. (|56l ) becomes 

1 /s = 0. (57) 



2^ 

dx 2 



a 2 2 
— T — X 

a 



Comparing Eqs. (|54]) and (0), one sees that a per4 plays the role of an effective background 
for the perturbations. An observer on the bulk brane sees a scale factor a B = D 1 / 2 (Eq. ©) 
but the fluctuations evolve, according to Eq. (|54"D, as if the scale factor were a pert . Hence, 
the shape of the fluctuation spectrum depends on a pert , not a B = D 1 ^ 2 ; but the physical 
wavelength is determined by a B = D 1 ^ 2 , not a per t- This is an important subtlety in our 
calculation. Let us now discuss the Hubble horizon for the perturbations. Recall that in 
usual 4d cosmology (see Eq. (f)7])), we have 



x = k(-r) = • a ■ (-t) = k phys a ■ (-r) ~ k phys H 1 , 



(58) 



where H^ 1 = a 2 /a is the Hubble radius as derived from the scale-factor a. By definition, 
a mode is said to be outside the Hubble horizon when its wavelength is larger than the 
Hubble radius. From Eq. (|58|), we see that this occurs when x < 1. Therefore, a mode with 
amplitude ft crosses outside the horizon when x ~ 0(1). Similarly, in our scenario we can 
write 

x = k{— t) = kphygDj (— r) = kphygHp^, (59) 

where fcphys = k/D^ 2 (since a# = .D,!/ 2 relates comoving scales to physical length scales on 
the bulk brane). The role of the Hubble radius is replaced by 



2{E-V{Y)) 
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which is to be thought of as an effective Hubble radius for the perturbations. So, as suggested 
above, the length scale at which amplitudes freeze depends on ag (rather than a pert ), but 
the amplitude itself, as derived from Eq. fl54"|), depends on a pert . The feature of two different 
scale factors is a novel aspect of the ekpyrotic scenario. 

By the time the bulk brane collides with the visible brane, modes are frozen on all 
scales less than the value of H~ ert when the bulk brane leaves the hidden brane. Comparing 
Eqs. (HH) and (|60|) , we see that this initial value of H~} rt is of the order of the particle 



horizon dnoR at the moment of collision. Recall from Section [IV D| that dnoR is required 
to be exponentially larger than the Hubble radius at collision, H' 1 , in order to solve the 
homogeneity problem. Hence, we see that modes are frozen on scales exponentially larger 
than if" 1 , thereby solving the inhomogeneity problem. 

(a) 
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FIG. 4. Sketch comparing the generation of a super-horizon spectrum of perturbations in (a) 
inflationary cosmology versus (b) the ekpyrotic universe. During inflation, the Hubble radius is 
nearly fixed and the fluctuation wavelength grows exponentially fast, causing modes to be stretched 
outside the horizon. In the ekpyrotic scenario, modes correspond to ripples on the moving bulk 
brane. The perturbations have nearly constant wavelength but the effective Hubble radius shrinks, 
once again causing modes to cross outside the horizon. 

The comparison to inflationary cosmology is made in Figure |j. The salient feature of 
both models is that perturbation modes inside the Hubble horizon escape outside in the 
early universe and re-enter much later. However, the behavior of the scale factor and the 
Hubble horizon are quite different. In inflation, the wavelengths are stretched superluminally 
while the horizon is nearly constant. In the ekpyrotic scenario, the wavelengths are nearly 
constant while the horizon shrinks. 

It remains to show that we can obtain a spectrum which is scale-invariant. Writing the 
equation for the perturbations in the form of Eq. (|54]) is useful since one can read off from 
it the spectral slope of the power spectrum. It is determined by the value of (d pert / 'a peTt )T 2 . 
In particular, one obtains a scale-invariant spectrum if (d per t/a per t)T 2 = 2 when the modes 
observed on the CMB cross outside the horizon. (Note that in usual 4d cosmology, this 
is achieved for an expanding de Sitter universe with a oc — r _1 or a contracting matter- 
dominated universe with a oc r 2 .) Combining Eqs. (52) and (|55|) , we find 



tipert 2 r->-3 ( ^Vq ^ d 2 Vt 

-r — D n a—— — D, 



2sx - * D(Y')dY' 



" 



2(E-V(Y')) 



(61) 



a per t \ dYo dY§ 

The spectrum will be scale-invariant if the right hand side of Eq. (|B"ID equals 2 when the 
modes of interest cross outside the horizon. 

As a simple example, consider the case where |V(y)| <C \E\. Eq. (|6"T| ) then generically 
gives (d pert /a pert )r 2 1. This leads to a density spectrum of the form |d> fc | ~ k, which 
is thus unacceptably blue. (A similar calculation was repeated for other set-ups such as 
Randall- Sundrum and the solutions presented in Ref. It was found that none of these 
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solutions predict a scale-invariant spectrum of fluctuations when V is turned off.) It is 
therefore crucial to add a potential in order to obtain a scale-invariant spectrum of density 
fluctuations. 



2. A successful example: The exponential potential 

We can add a potential V(Y) of the form that might result from the exchange of wrapped 
M2-branes. We would like to think of V as the potential derived from the superpotential W 
for the modulus Y in the 4d low energy theory. Typically, superpotentials for such moduli 
are of exponential form, for example, 

W ~ e~ cY , (62) 

where c is a positive parameter with dimension of mass. The corresponding potential is 
constructed from W and the Kahler potential K according to the usual prescription 



V = e K ' M ; 



K ij DiW D~W — tptWW 
M pI 



(63) 



where D$ = d/d<f) 1 + Ki/M% is the Kahler covariant derivative, Ki = dK/dfa, Kij = 
d 2 K / d(j? d<pi and a sum over each superfield fa is implicit. Eqs. ( [52] ) and imply that V 
decays exponentially with Y. For the purpose of this paper, we shall not worry about the 
exact form of the superpotential in heterotic M-theory. Rather, it will suffice to perform the 
calculation using a simple exponential potential, namely 

V(Y) = -ve- maY , (64) 

where v and m are positive, dimensionless constants. In this paper, it is convenient to 
parameterize the exponent in terms of a. Note that, in the case where the potential is 
generated by the exchange of wrapped M2-branes, the parameter m is of the form m = 
cT^u/a, where c is a constant, T3 is the tension of the M2-brane, and v is the volume of the 
curve on which it is wrapped. The potential defined in Eq. ( p^ ) is shown in Figure 0. The 
perturbation modes of interest are those which are within the current Hubble horizon. As 
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the wavelengths corresponding to those modes passed outside the effective Hubble horizon on 
the moving bulk-brane, the amplitudes became fixed. Scale invariance will require mD ^> 1 
during this period. (In section |V B 3| , we shall generalize this condition for an arbitrary 
potential V.) 

We have already seen at the end of Section |V B 1| that, if the potential V is negligible 
compared to E, the spectrum of fluctuations is not scale- invariant. Hence, we consider the 
limit where \E\ -C |Vq|. This condition, as seen from the equation of motion for Yq, Eq. ([51~1), 
is satisfied if Yq = initially, or, equivalently, if the bulk brane begins nearly at rest. For 
the brane to be nearly at rest, one must have \E\ « \Vq\ initially. As the brane traverses the 
fifth dimension, \V\ increases exponentially, whereas E is constant. Hence, the condition 
\E\ <C | Vo| is automatically satisfied. The bulk brane beginning nearly at rest is precisely 
what we expect for a nearly BPS initial state. 

Applying the condition \E\ -C |Vq|, Eq. (ET) reduces to 



2 1 

~ 2v 



[ Y ° D(Y')e maY ' /2 dY' 
Jo 



2D 2 



o 



m 2 a 2 ve maY o \ mD 



(65) 



where we have neglected the endpoint contribution at Y — 0. On the other hand, Eq. (| 
gives 

tipen = m 2 a 2 ve- maYo ( 1 \ 
apert E>1 V mD J ' 

Combining the above two expressions yields 

^r 2 = 2(l + -*—) ( 1 - —) . (67) 
a pert \ mD J \ itiDqJ 

The right hand side of Eq. ( |6^) is approximately equal to 2 in the limit of large itlDq. 

Hence, the exponential potential of Eq. ( j64| ) results in a nearly scale-invariant spectrum 

of perturbations provided that \E\ |Vo| and mD ^> 1 are satisfied when modes pass 

outside the effective Hubble horizon. Note that it would be exceedingly difficult to maintain 

(tipert I 'apert)T 2 close to almost any other value than 2. It is indeed fortunate that scale- 

invariance is the desired result because obtaining a different spectral index from the ekpyrotic 

scenario would be highly problematic. 
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We next compute the perturbation amplitude, by using Eq. ([54]) to calculate |Alfc|. As 
shown above, in order for the spectrum to be scale- invariant, the conditions \E\ -C |Vq| and 
itlDq ^> 1 must be satisfied when wavelengths pass outside the horizon. These conditions 
can be relaxed once the mode is well outside the horizon. For example, we will assume no 
restrictions on mC, the value of mD at y = 0. 

In the limit that mDo 3> 1 when the relevant modes cross outside the horizon, Eq. (|54"1) 
reduces to 



X 



dx 2 



x 2 



H = 0, (68) 



with solution 



f k = x 1 ' 2 (C^k) J 3/2 (x) + C 2 (k)J- 3/2 (x)) , (69) 

where J±s/ 2 are Bessel functions. The coefficients C\{k) and C 2 (k) are fixed by requiring 
that modes well- within the horizon (i.e., x ^> 1) be Minkowskian vacuum fluctuations, that 



is 



h = e' lkT for x » 1. (70) 

/6kpM$B 



(Note that the factor of y6/3Mf.B in Eq. (|70|) arises when we change variables from Y to a 
canonically- normalized scalar field.) Using this initial condition, we find the following rms 
amplitude for modes outside the horizon (with s<l) 

a f — k 3/2 fk -i /y.N 

Substituting Eq. (|65D and using fk = D 5Yk, we find 



Finally, we define the time-delay Ar(k) by 

AY k 



\Ar(k)\ 



ma \/ve~ maYo 

An = V n2 • (72) 

2(2tt) 3 / 2 v / 3 / 3M 5 3 5 d o 



Yn 



167r 3 / 2 y / 3^Mp ^ mD 
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where we have used the equation of motion for Y , Eq. (0). Note that the time-dependence 
of Ar(k) is mild, a necessary condition for the validity of the time-delay formalism. The 
factor of mDo = mD(Yo(r)) is to be evaluated at time r when a given mode crosses outside 
the horizon during the motion of the bulk brane. Let Dk denote the value of Do at horizon 
crossing for mode k. Since horizon crossing occurs when x = 1, or, equivalently, when 
(— r) = k^ 1 , Eq. (|65| ) gives 



m \ 2k 

Substituting Eqs. (^) and @ into Eq. @, we find 



oft«VS / 2 



47r 3 / 2 v /3/3M 5 3 J B(mC + 2) VmA 

This expression for \Sk\ increases gradually with increasing k, corresponding to a spectrum 
tilted slightly towards the blue. The blue tilt is due to the fact that, in this example, D is 
decreasing as the brane moves. That is, the spectral index, 

d log \5k | 2 _ 4 

™ s = 1 + "Ji — i - ~ 1 + — FT' 76 
a log k mlJi 

exceeds unity. The current CMB data constrains the spectral index to lie in the range about 
0.8 < n s < 1.2. Therefore, for our results to be consistent with experiments, we must have 

mD k > 20, (77) 

a constraint that is easily satisfied. 

3. General potential 
As a second example, consider the power-law potential 

V(Y) = -vD(Y) q = -v{aY + C) q , (78) 
where v > and q < are constants. In this case, Eq. (^T|) gives 
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Q'pert ^_2 ^ 2 



1-4 

9 



(79) 



for |g| ^> 1. Hence, a power- law potential can also lead to a nearly scale-invariant spectrum 
provided that its exponent is sufficiently large. (The smaller is the value of |g|, the bluer is 
the spectrum.) 

We can straightforwardly extend our analysis to an arbitrary potential V{Y). Let us 
suppose that V(Y) satisfies 



>a\V(Y)\ 



D{Y) 

(For the exponential potential, V(Y) - 
Then, Eq. (|H]) reduces to 



d 2 V 



dY 2 





dV 


> a 


dY 



-ve 



(80) 



maY these two conditions amount to mD 3> 1. 



a 



pert 



(I 



pert 



'VV' r 
y/2 



(81) 



Hence, the conditions for scale invariance are Eqs. (|80|) as well as 

VV" 

w 1. 

y/2 



(82) 



C. Gravitational waves from colliding branes 

In inflationary cosmology, the analysis of tensor (gravitational wave) perturbations fol- 
lows closely the analysis of scalar (energy density) perturbations.0 Metric fluctuations can 
be divided into two polarizations, each of which acts like a massless scalar field evolving in 
the same cosmic background as the inflaton. Hence, it is not surprising that the spectrum of 
tensor fluctuations has nearly the same scale-invariant spectral shape as the scalar spectrum. 

In the ekpyrotic scenario, the relationship between scalar and tensor perturbations is less 
direct. The excitations that produce scalar perturbations are ripples on the moving brane, 
which are directly dependent on the rate at which the brane traverses the fifth dimension 
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and the potential that drives it. The tensor fluctuations, on the other hand, are excitations 
of the gravitational field, which lives in the bulk. (The moving brane itself does not support 
tensor fluctuations.) The net result is a different effective scale factor in the fluctuation 
equation of motion for the tensor modes than for the scalar modes. 

We shall briefly outline the derivation here, with more details to follow in our more formal 
paper on perturbations.^ If is the unperturbed, homogeneous metric (see Eq. (|ll]) with 
A and N functions of time), the perturbed 5d metric can be written as 

9^ = 9^u + A 2 (t)D(y, t) V(x, t), (83) 

where /i, v = 0, . . . 3. Note that since we shall work at the level of the 4d effective theory, 

we can treat the tensor perturbations h^ u as functions of x and t only. We are interested in 

the tensor perturbations which satisfy the conditions: h 0fJi = 0, hj = 0, and d l hij = 0. The 

perturbed 5d Einstein action to quadratic order is 

y\/f3 r y\/f2 t 

Sj luc t = d'x^-gR d 4 xa 2 (h^h» - d t h^h») (84) 

where the second expression is obtained by integrating over y. The tensor action is analogous 
to the scalar action given in Eq. flo*5|). From the action, we can derive the tensor analogue 
of the scalar fluctuation equation of motion, Eq. (|54|) 



d 2 fT 
2 a Jk 



dx 1 

where 



a 2 2 
-T — X 



ft = 0, (85) 



d 3 k 

K = I ^ ( "« ,,>AT) (86) 



and 



H=^H- ( 87 ) 

The critical difference between this tensor equation and the scalar fluctuation equation, 
Eq. fl54"D, is that the effective scale factor a pert in Eq. fl54]) has been replaced by a. 
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We introduced a potential to insure that a pert led to a nearly scale-invariant spectrum, 
iflpertl 'dpert)' 1 ' 2 ~ 2. However, a(r) in the tensor equation is approximately constant (recall 
that a = (BI^Ms) 1 / 2 + O(0/ a)). Consequently, the root mean square tensor fluctuation 
amplitude, 

k 

' A ^' ^ww~ww (88) 

is not scale- invariant. Rather, the tensor spectrum is tilted strongly to the blue. Fitting 
the mean square amplitude to a scale- free form, ~ k UT , where tit is the conventional tensor 
spectral index, the spectrum above corresponds to nx = 2, compared to the inflationary 
prediction, tit < 0. 

The tensor spectrum is a prediction that clearly distinguishes the ekpyrotic scenario 
from inflationary cosmology. In both cases, for the same Hubble parameter at reheating (for 
inflation) or collision, H c , the mode with wavelength of order if" 1 has similar amplitude, 
H c /M p i, where M pi is the 4d Planck mass. The wavelength of this mode today is roughly 60 
e-folds smaller than the current Hubble radius. Hence, if we extrapolate from this wavelength 
to one comparable to the present Hubble radius, Hq 1 , the inflationary prediction is that the 
amplitude is nearly the same (since the spectrum is nearly scale-invariant), whereas the blue 
spectrum computed above predicts that the amplitude is exponentially small. Hence, the 
search for a gravitational wave signal using the CMB polarization on horizon scales becomes 
a key test for our proposal. Future gravitational wave detectors, beyond the presently 
planned LISA and LIGO projects, may also someday detect the stochastic background of 
gravitational waves as well. Observing a nearly scale invariant primordial gravitational wave 
background falsifies the ekpyrotic scenario and is consistent with inflation. 



VI. CONCLUSIONS 
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A. Recapitulation 



Conceptually, the ekpyrotic scenario appears to be simple: a bulk brane strikes our visible 
brane and a hot big bang universe is born. In actuality, to accomplish the transformation 
from a cold, nearly BPS state into an expanding, hot universe with nearly scale-invariant 
density perturbations without invoking inflation requires a series of seemingly incongruous 
conditions. Remarkably, these conditions can be satisfied simultaneously in heterotic M- 
theory. 

First, the gravitational backreaction due to the kinetic energy of the bulk brane must 
trigger cosmic expansion. In 4d gravity, kinetic energy usually causes cosmic deceleration 
and, if the initial state is static, it triggers contraction. Second, the total energy of the brane 
has to grow by drawing energy from the gravitational field since, otherwise, the total bulk 
brane energy before and after collision is zero and there is no radiation. In 4d gravity, this 
blue shift effect occurs if the universe is contracting, but here it occurs even though the scale 
factor on the visible brane is expanding. Third, the scalar (energy density) fluctuations must 
be nearly scale- invariant. Although scale-invariance is ordinarily associated with inflation, 
here we have shown that a scale-invariant spectrum results even though the universe is 
quasi-static. The only requirement is a bulk brane potential whose magnitude increases 
by an exponential factor as the brane traverses the bulk. Potentials of this type occur in 
string theory. (Our analysis suggests that scale-invariance is especially favored in the 5d 
theory. It occurs for rather simple, physically-motivated potentials, whereas more general 
spectral shapes are difficult to obtain.) Fourth, the tensor fluctuation spectrum is not scale- 
invariant but, rather, strongly tilted towards the blue, providing an observational signature 
that distinguishes the ekpyrotic scenario from inflationary cosmology. 

In the context of 4d gravity, some of these features suggest slow expansion, others super- 
luminal expansion, yet others contraction. How do we obtain all of these features simulta- 
neously? All of this is possible in our 5d theory because the role of gravity in the equations 
of motion for the expansion, the brane motion, and the scalar and tensor fluctuations is 
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played by different combinations of moduli fields for each equation. That is, each equation 
is similar to the corresponding equation for a scalar field in a 4d FRW background except 
that the scale factor a(t) is replaced by some function of the moduli fields that differs for 
each physical quantity. Some combinations increase with time, mimicking an expanding uni- 
verse, and others decrease, mimicking a contracting universe. The remarkable result is that 
just the right combinations of moduli occur in heterotic M-theory to produce the behavior 
required for a viable cosmological scenario. 

An observer at any surface of fixed y has a scale factor equal to AD 1 ^ 2 (y, r) (see Eq. (|TTD ), 
which is expanding as the bulk brane collides with the visible brane. In Eq. fl23l) , which 
describes the time variation of the total energy of the bulk brane, the role of the scale factor 
is played by a, which is contracting. The contraction produces the blue shift or increase in 
the total energy so that, upon collision, there is excess kinetic energy that can be converted 
to radiation. The scalar fluctuations are ripples in the bulk brane surface that evolve as if 
the scale factor were a pert in Eq. fl54"|), corresponding to a contracting effective Hubble radius. 
We have identified simple criteria for the potential which result in a scale- invariant spectrum. 
The tensor fluctuation spectrum naturally differs from the scalar spectrum because tensor 
fluctuations occur in the bulk volume rather than on the brane surface. The effective scale 
factor for the tensor fluctuations is a, rather than a pert . The differences account for the fact 
that the scalar spectrum is scale-invariant, whereas the tensor spectrum is tilted strongly to 
the blue. 

A useful mnemonic for recalling the difference between the scalar and tensor fluctuation 
spectra in our scenario is to consider the equivalent relations for inflationary cosmology, but 
with the inflaton scalar field replaced by Y. For the scalar fluctuations, the amplitude is 




(89) 



and the tensor fluctuation amplitude is 




(90) 
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where M v \ is the 4d Planck mass, Ay is the fluctuation amplitude for Y, and the subscript 
means that the expressions are to be evaluated when the wavenumber of a given mode 
is equal to the inverse Hubble radius H as it passes beyond the horizon. For inflation, 
H, AY and Y are nearly constant, so 5 s and 5t are both nearly scale-invariant. For the 
ekpyrotic scenario, Ay and Y are both strongly time-varying. However, for an exponential 
potential V(Y) as is naturally generated by non-perturbative exchange of M2-branes, the 
time-variation in the ratio AY/Y in the expression for 5s nearly cancels. Consequently, the 
ratio is nearly constant and the resulting spectrum is nearly scale-invariant. However, 5t 
involves only H, which is increasing with time as smaller and smaller wavelength modes 
pass beyond the horizon. This accounts for the fact that the spectrum is blue. 

As a specific fully-worked example, consider the case of an exponential bulk brane po- 
tential, V(Y) = —vexp(-maY), as discussed in Section [VB 2| . We have computed the 
ekpyrotic temperature at the beginning of the hot big bang phase (Eq. (|30|)), 

_T_ tf/\2v) 1 ^ ( Mg \ 1/2 f(3\ 1/2 (mC + 2) 1 / 2 

M pl ~ (I 3 M 5 y/*(aR + Cy/* \M pl ) [a) m ( } 

In terms of this temperature, the scalar (energy density) fluctuation amplitude in Eq. ([75]) 
can be rewritten as 

\5 1 = m ^ haf ' 2 (°) 3/2 -A_ (Z\ 2 m) 

1 *' 367r 3 / 2 (mC + 2) 2 \(3J mD k \M pl J ' { ' 

A simple example which satisfies all constraints is a = 2OOOM5, (3 = M5, B = 10~ 4 , 
C = 1000, R = Ms" 1 , m = 0.1, and v = 10" 8 , all of which are plausible values. In 
this example, (the value of D at horizon crossing) is of order 10 3 . Then, we find that 
M5 ~ 10~ 2 M p ;; the ekpyrotic temperature is T ~ 10~ 8 M p ;; and the scalar perturbation 
amplitude is \5k\ ~ 10~ 5 . Note that the ekpyrotic temperature, the maximal temperature 
of the hot big bang phase, tends to be small compared to the Planck or unification scale. 
This is a characteristic feature of the model. With these parameters, the magnitude of 
the potential energy density at collision is (10 _6 M p ;) 4 . This corresponds to a characteristic 
energy scale for the potential of 10 13 GeV. Finally, note that these values are consistent 
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with Hofava-Witten phenomenology. For instance, the proper distance between the branes 
is Rproper ~ 10~ 5 M p ;. If we further assume that the characteristic length scale L C y of the 
Calabi-Yau three-fold is approximately 10 times smaller than R pr0 per (in order for the five- 
dimensional effective theory to be valid), then we get agreement with the values of R pr oper, 



Ley, and the lid Planck mass M\\ inferred by Witten in Ref. ITS] in matching the gauge 
and gravitational coupling constants. 

While we are pleased that the numerical constraints from cosmology and those from 
Hofava-Witten phenomenology can be simultaneously satisfied, we should emphasize that 
there is a lot of flexibility in terms of parameters. For instance, in the above example the 
ratio (3 /a was chosen to be of order 10~ 3 . However, one can easily make this ratio as large 
as 1/10 if one wishes. For example, choosing a = 200M 5 , j3 = 20M 5 , B = 10~ 2 , C = 100, 
R = M 5 -\ m = 1, and v = 10~ 10 results in T ~ 10~ 7 M ph \5 k \ ~ 10" 5 , M 5 ~ HT 2 M ph 
Rproper ~ 10" 5 M p z, and characteristic energy scale for the potential of 10 14 GeV. The value 
of Rproper could also take a significantly different value, if one wishes, and there would still be 
enough freedom to obtain reasonable ekpyrotic temperature and fluctuations. We can also 
imagine applying the same ideas in a different brane world context, such as AdS, and still 
obtaining a successful scenario from a cosmological point-of-view. The challenge, of course, 
is to figure out how to break supersymmetry, obtain a correct phenomenology and stabilize 
moduli. Here we have presumed that this challenge can be met, and have shown through 
examples how the brane world approach to particle phenomenology might be combined with 
new ideas in cosmology to obtain a successful picture of the early Universe. 

B. Colliding Branes and Inflationary Cosmology 

We have laid out a detailed cosmological scenario that offers a resolution of the flatness, 
horizon, and monopole problems and generates a nearly scale-invariant spectrum of energy 
density perturbations based on concepts that derive naturally from extra dimensions, branes, 
and heterotic M-theory. The key conceptual difference from inflation is how the universe 



56 



begins. In the usual approaches to inflationary cosmology, as in standard big bang cosmol- 
ogy, the universe begins with a cosmological singularity. The universe emerges in a high 
energy state with no particular symmetry and rapidly expanding. Superluminal expansion 
is invoked to smooth out and flatten the emerging state. The ekpyrotic scenario intro- 
duces a different philosophy in which the universe begins in a non-singular, infinite, empty, 
quasi-static state of high symmetry. Superluminal expansion is not needed because the BPS 
vacuum state is flat and smooth. Brane collision can account for the matter-radiation energy 
and primordial density perturbations. 

Let us briefly summarize how the ekpyrotic scenario addresses the various cosmological 
problems: 

• Causal Horizon Problem: In the ekpyrotic scenario, the local temperature and density 
are set by the collision of the visible brane and bulk brane, which acts as a non-local 
event that occurs nearly simultaneously over a region much larger than the Hubble 
horizon. 

• Flatness Problem: The universe is assumed to begin in a nearly BPS ground state. 
The BPS state corresponds to a spatially-flat geometry. The process of bulk brane 
formation/nucleation and propagation maintains flatness. (We do not demand that 
the initial state be globally BPS to resolve the horizon and flatness problems. It 
suffices that the universe be flat and homogeneous on scales up to the (causal) particle 
horizon, as should occur naturally beginning from more general initial conditions. In 
the ekpyrotic scenario, because the bulk brane motion is extremely slow, the particle 
horizon at collision is exponentially large compared to the Hubble horizon, where the 
latter is set by the radiation temperature after collision.) 

• Monopole Problem: The hot big bang epoch commences when the bulk brane col- 
lides with the visible brane and heats the universe to a finite temperature. Provided 
the temperature is less than the monopole mass, the monopole abundance will be 
negligible. 
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• Inhomogeneity Problem: Quantum fluctuations generate ripples in the bulk brane 
as it traverses the bulk. Due to the ripples, collision and thermalization occur at 
varying times across the visible brane, resulting in fluctuations in energy density and 
gravitational waves. 

Both the ekpyrotic scenario and inflationary cosmology have the feature that the causal 
horizon is exponentially greater than the Hubble horizon. In inflation, superluminal expan- 
sion rapidly stretches the causal horizon while the Hubble horizon is nearly fixed. In the 
colliding brane picture, the collision of the bulk brane acts as a non-local interaction that 
causally links regions separated by much more than a Hubble distance. 

Both the inflationary and ekpyrotic scenarios produce a nearly scale-invariant spectrum 
of energy density perturbations from quantum fluctuations. For inflation, quantum fluc- 
tuations are stretched beyond the Hubble horizon as the universe expands superluminally. 
For the ekpyrotic universe, the Hubble horizon is shrinking compared to the quantum fluc- 
tuations as the universe contracts very slowly. The equations describing the evolution of 
perturbations are (nearly) equivalent in the two cases (see discussion of Eq. (0)), even 
though, one describes an expanding de Sitter phase and the other a contracting pressureless 
phase. The similar equations account for why both lead to scale-invariant spectra for density 
perturbations even though the mechanisms are different. 

From the point-of-view of an observer on the stationary orbifold planes, the universe is 
expanding as the branes collide. The bulk brane is what causes their expansion, a grav- 
itational backreaction effect due to its motion. The expansion is very slow as the brane 
moves across the fifth dimension, but assumes the usual big bang rate after collision and 
thermalization. 

One might hope that the ekpyrotic scenario avoids the tuning problems required in stan- 
dard inflation in order to obtain an acceptable perturbation spectrum. Thus far, the situa- 
tion is unclear. We found that we had to introduce a flat potential for the bulk brane that 
is roughly similar to the flat inflaton potential used in standard inflation. The form is also 
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qualitatively consistent with non-perturbative potentials that arise in M-theory. Perhaps the 
potential parameters needed for our scenario will be shown to arise naturally. However, it 
should also be noted that the reasons for introducing the potential in the ekpyrotic scenario 
are different from the case of inflation. In our case, the need for a flat potential is linked 
to the precise form of the background static BPS solution of Lukas, Ovrut and Waldram 
used in this paper. Perhaps there exist other initial conditions which avoid the need for flat 
potentials altogether. 

Although inflationary cosmology and the colliding brane picture both produce a nearly 
scale-invariant spectrum of perturbations, the deviation from scale-invariance differs due 
to the fact that the background felt by the perturbations is expanding in one scenario 
and contracting in the other. In standard inflationary cosmology, the spectrum of scalar 
(density) and tensor (gravitational wave) perturbations is typically red (amplitude decreases 
as wavelength decreases) .0 The amplitude of a given mode is proportional to the Hubble 
parameter when the wavelength is stretched beyond the horizon. The Hubble parameter 
decreases (slowly) in an expanding, inflating universe. Since smaller wavelength modes 
stretch beyond the horizon at later times when the Hubble parameter is smaller, they have 
a smaller amplitude, resulting in a red spectrum. The degree of redness is expressed in 
terms of a "spectral index, '@ tis,t for scalar and tensor perturbations respectively, where 
rig — 1=71,^ = is defined as precise scale-invariance, and rig — 1 < and < correspond 
to red spectra. In the examples of the ekpyrotic scenario discussed here, the apparent Hubble 
radius for an observer on the bulk brane is shrinking. Consequently, the corresponding 
spectra are blue. By introducing a potential for the bulk brane (dependent, say, on its 
position Y), the density perturbation spectrum can be made nearly scale-invariant, slightly 
blue (ns — 1 > 0) in our examples. On the other hand, the gravitational wave spectrum is 
unaffected by the potential and is strongly blue riy ~ 2. 

For energy density perturbations, there are exceptional cases where inflation can give a 
blue spectrum.il The blue spectrum arises because the density perturbation amplitude is 
not only proportional to the Hubble parameter, but also inversely proportional to the kinetic 
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energy of the inflaton. As inflation proceeds, the Hubble parameter decreases and, in most 
models, the inflaton kinetic energy increases; so both effects tend to make the spectrum red. 
But models can be rigged where the Hubble parameter decreases, as usual, but the inflaton 
kinetic energy decreases more rapidly In that case, the spectrum is blue. Similarly, it is 
possible to get a red spectrum in the ekpyrotic model, for example, if the bulk brane moves 
in the direction of increasing warp factor. Hence, observing a red or blue density spectrum 
is not a decisive test for distinguishing the two scenarios. 

However, the gravitational wave spectrum for inflation is always red — the amplitude 
depends only on the Hubble parameter — and so observing a strongly blue gravitational 
wave spectrum, as predicted by the ekpyrotic scenario, is a key test. The cosmic microwave 
background polarization is one method of detecting the presence of primordial gravitational 
waves with wavelengths comparable to the Hubble horizon today. For the slightly red spec- 
trum of inflationary cosmology, the gravitational wave amplitude may be large enough to be 
detected. However, for a strongly blue spectrum, the signal at large length scales is exponen- 
tially small and undetectable. Hence, the detection of gravitational waves in the microwave 
background polarization would falsify our scenario. Currently planned gravitational wave 
detectors, such as LIGO and LISA, are not sensitive enough to detect the gravitational waves 
from inflation or from our scenario. However, future detectors may discover the stochastic 
background and determine the spectral slope. 

Certain aspects of our scenario are reminiscent of the novel, string-inspired, pre-big bang 
models introduced by Veneziano and Gasperiniilllll Both assume the universe begins in a 
flat, empty state.@ In both models, the gravitational wave perturbation spectrum is blue.ll 
However, the structure, ingredients, dynamics and predictions of the two models are very 
different. The pre-big bang scenario does not entail extra dimensions or branes in a direct 
way. The pre-big bang begins with a semi-infinite period of contraction which is superluminal 
(deflation). This period ends in a global singularity in which the Hubble constant becomes 
infinite. Matching across this singularity is the biggest challenge facing the pre-big bang 
model. If it is possible at all, it can only be at the string scale where non-perturbative stringy 
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effects are large and difficult to compute. In our model, matching to the universe after the 
bulk-boundary collision is complex, and will require calculations in the five dimensional 
theory which we have not studied here. But the collision event itself is nonsingular and 
it is an important advantage of our scenario that it only involves physics taking place at 
low energies, which is in principle describable using the effective low energy Lagrangian. A 
second difference is that all expansion and contraction are subluminal in our model. Finally, 
in the simplest renditions of pre-big bang, the fluctuation spectrum is so strongly tilted to 
the blue relative to a scale invariant spectrum that current observational bounds on the 
microwave background anisotropy are violated. The ekpyrotic scenario obtains a spectrum 
that is consistent with current observations. 

To summarize, we have presented a novel scenario for the beginning of the hot big bang 
universe, within a framework consistent with string theory and supergravity. The universe 
begins in the simplest state possible, one which is cold, nearly BPS, and nearly vacuous. At 
some time, a bulk brane exists or is nucleated in the vicinity of the hidden brane (through 
a small instanton phase transition), and begins to move towards the visible brane. The 
bulk brane eventually collides with the visible brane and is absorbed in a small instanton 
phase transition. This transition may change the gauge group on the visible brane to the 
standard model gauge group, as well as create three families of light quarks and leptons. At 
the moment of collision, a fraction of the kinetic energy of the brane is converted to thermal 
excitations of the light degrees of freedom on the visible brane, causing the universe to 
enter an FRW radiation-dominated phase. Furthermore, ripples on the bulk brane imprint 
a spectrum of energy density fluctuations consistent with current observations and which 
provides the seeds for structure formation. While parts of our scenario remain speculative 
at present (such as the dynamics of the small instanton phase transition), it is our hope 
that advances in heterotic M-theory will eventually allow us to solidify the components of 
our cosmological model. For the moment, we consider our scenario as a first step towards a 
new, testable model for the early universe consistent with current cosmological observations 
and fully-motivated by string theory. 
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